V 


A~ 


AN  INVESTIGATION  OF  THE  EFFECT  OF 
AUTO-THROTTLE  DEVICES  ON  AIRCRAFT  CONTROL 
IN  THE  CARRIER  LANDING  APPROACH 


*   *   *   *   * 


Gerald  R.  Bell 


AN  INVESTIGATION  OF  THE  EFFECT  OF 

AUTO-THROTTLE  DEVICES  ON  AIRCRAFT  CONTROL 

IN  THE  CARRIER  LANDING  APPROACH 


by 

Gerald  R.  Bell 

// 


Lieutenant  Commander,  United  States  Navy 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 
AERONAUTICAL  ENGINEERING 


United  States  Naval  Postgraduate  School 
Monterey,  California 


1963 


DUDLEY  KNOX  UBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY  CA  93943-5101 


if 


NAVAL  I  E SCHOOL 

AN  INVESTIGATION  OF  THE  EFFECT  OF 

AUTO -THROTTLE  DEVICES  ON  AIRCRAFT  CONTROL 

IN  THE  CARRIER  LANDING  APPROACH 

by 

Gerald  R.  Bell 


This  work  is  accepted  as  fulfilling 
the  thesis  requirements  for  the  degree  of 
MASTER  OF  SCIENCE 
IN 
AERONAUTICAL  ENGINEERING 
from  the 
United  States  Naval  Postgraduate  School 


ABSTRACT 

An  automatic -throttle  compensation  system  has  the  capability 
of  removing  the  problem  of  aircraft  speed  instability  in  the  carrier 
landing  approach.  Various  systems  have  been  proposed  and  tested,  each 
differing  in  input  sensed  variables.    Two  representative  systems  are  in- 
vestigated analytically  and  by  the  use  of  digital  programmed  Nyquist 
plots  and  analog  simulations.    An  attempt  is  made  to  determine  optimum 
gain  constants  by  using  various  forcing  functions  in  the  analog  simulation. 
Comparisons  of  the  responses  of  the  two  systems  are  made  using  time 
history  analog  records.    A  digital  computer  stability  program e  applicable 
to  any  aircraft,  is  included. 
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AN  INVESTIGATION  OF  THE  EFFECT  OF  AUTO-THROTTLE  DEVICES 
ON  AIRCRAFT  CONTROL  IN  THE  CARRIER  LANDING  APPROACH 


1 .         Introduction 

The  carrier  landing  approach  speed  of  current  swept  wing  aircraft 
is  influenced  by  many  factors .    One  of  the  most  important  of  these  factors 
is  often  referred  to  as  speed  instability.    Speed  instability  can  be  ex- 
plained with  the  help  of  Fig.  1 .    It  will  be  noted  that  for  Aircraft  "A"  ,  a 
slight  deviation  in  airspeed  from  the  normal  approach  speed  makes  very 
little  difference  in  thrust  required.    This  curve  is  typical  of  the  AF-1E 
type  aircraft.    On  the  other  hand,  a  deviation  towards  slower  speed  for 
Aircraft  "B"  causes  a  large  increase  in  thrust  required.,    This  increase , 
if  uncorrected,,  causes  a  loss  of  altitude„    Thus,  in  effect,  the  normal 
approach  of  Aircraft  "B"  is  made  on  the  "backside"  of  the  thrust  versus 
airspeed  curve      This  curve  is  typical  of  the  F-8  aircraft „    Essentially 
this  steep  backside  of  the  curve  is  caused  by  a  high  value  of  Cn      , 
drag  increase  due  to  angle  of  attack  (or  lift) . 

Still  another  factor  contributing  to  airspeed  instability  is  slow 
engine  response.    With  a  comparatively  long  time  lag  in  engine  response, 
any  hesitation  on  the  part  of  the  pilot  in  correcting  an  incipient  airspeed 
loss  leads  to  wide  variations  in  altitude  or  airspeed  before  effective 
thrust  can  be  realized. 

With  an  aircraft  that  displays  speed  instability  tendencies ,  either 
a  faster  approach  speed  or  an  automatic  compensating  device  is  required 
to  prevent  undesired  altitude  loss  due  to  airspeed  deviations.    Faster 


approach  speeds  are  not  desirable  for  carrier  approaches.    It  will  be  the 
purpose  of  this  paper  to  investigate  some  methods  of  applying  automatic 
power  compensation  to  counteract  airspeed  instability .    A  swept  wing 
fighter,,  patterned  after  the  F-8  airplane  is  used  as  the  vehicle  of  the 
investigation.    Both  digital  and  analog  computers  are  exploited  in  simu- 
lating auto-throttle  and  anti-drag  system  response  to  various  inputs  „ 
Real-time ,  pilot  controlled,  analog  simulated  mirror  approaches  are  also 
used  in  qualitatively  evaluating  each  system  , 

This  investigation  is  the  first  of  a  continuing  series  to  be  conducted 
at  the  U.S.  Naval  Postgraduate  School  to  evaluate  systems  of  automatic 
airplane  control . 


2.         Discussion 

The  aircraft-carrier  mirror-landing  approach  is  made  at  a  constant 
speed  and  a  constant  angle  of  glide  slope „    This  type  of  approach  allows 
precision  control  of  the  landing,  a  factor  of  the  utmost  importance  in 
carrier  operations.    However ,  the  approach  is  usually  made  at  a  minimum 
airspeed  that  is  dictated  by  performance  and/or  aircraft  flying  qualities. 
Thus,  the  airspeed  must  be  monitored  closely  while  making  flight  path 
angle  adjustments ,    A  simplified  equation  (given  in  Ref.  1)  shows  that 
flight  path  angle,  ^  ,  is  approximated  by  the  formula: 

EC.  (i)  *=L      &+JHJ         ■        -.,    ^ 

Flight  path  angle,  then,  is  dependent  on  lift/drag  (L/D)  ratio,  thrust/ 
weight  ratio ,  and  rate  of  speed  change .    In  the  mirror  approach  the  pilot 
can  adjust  V  and  airspeed  by  manipulating  the  throttle  and  elevator.    The 
manner  in  which  each  of  these  controls  can  be  used,  however,  is  depend- 
ent upon  the  L/D  ratio „    For  aircraft  with  a  high  L/D  ratio,  the  elevator  is 
used  to  control  ^  ,  while  thrust  is  used  to  maintain  constant  airspeed. 
With  a  L/D  ratio  in  the  vicinity  of  4  -  5  the  rate  of  speed  change  is  large 
during  maneuvers.    If,  in  addition,  the  L/D  ratio  decreases  with  increasing 
CL  in  the  approach  speed  range  then  speed  instability  is  accentuated, 

A  variation  of  lift/drag  ratio  versus  lift  coefficient  for  the  F-8  air- 
plane is  shown  as  Fig,  2.    It  can  be  seen  that  at  the  approach  speeds 
used,  (  1.1  -  lo2  VgL  ) ,  L/D  decreases  rapidly  with  increase  in  CL „ 
As  reported  in  Ref.  1,  pilots  stated  that  in  order  to  change  *£  for  an 


airplane  of  this  class,  increased  reliance  on  thrust  rather  than  elevator 
was  necessary.    In  the  case  of  the  F-8  airplane,  however,  height  response 
to  throttle  input  is  slow  due  to  small  thrust  line  angle  of  attack  and 
negligible  trim  changes  due  to  power.    The  F-8  pilot,  therefore,  is  forced 
into  using  the  stick  as  the  primary  flight  path  control „    Then,  since  the 
airspeed  mode  of  oscillation  (the  phugoid)  is     lightly  damped,  the  aircraft 
will  hunt  about  a  new  equilibrium  airspeed  and  glide  slope  angle  in  res- 
ponse to  an  elevator  deflection.    These  speed  changes  in  the  5  to  10 
seconds  following  elevator  input  are  disturbing  to  the  pilot. 

Another  way  of  interpreting  airspeed  instability  is  referred  to  in 
Ref.  2.    A  curve  similar  to  Fig0  1  can  be  drawn  for  thrust  required  versus 
dynamic  pressure,  q.    A  definition  of  a  stable  speed  regime  can  then  be 
postulated  as  that  portion  of  the  curve  where 

Ineq.  (2)  d  T^t^">        s,  Q 

To  express  this  inequality  in  terms  of  the  lift  and  drag  coefficients, 
we  use  the  fundamental  equilibrium  equations  for  level  flight: 

ClS<j.    =Lift   =  w 

therefore  T*.  -  W  ^& 


and 


since  W  and  S  are  constant. 


Therefore  d  TVgq.  £  J  (  ^  J 

which  result  is  substituted  into  Ineq.  (2)  .  S        (  Cl  /    s. 

Since  S  is  greater  than  O,  this  implies  that 
Taking  the  differentials  of  numerator  and  denominator: 


or 


Ineq.  (3);  C©'_   */£& '     >   Q 

Thus,  satisfaction  of  Ineq .  (3)  implies  speed  stability. 

The  C    °  term  used  above  is  composed  of  C     +  CAC/  where  C.„  is 
D  D         AS'  AS 


equivalent  propulsion  system  drag.    C._,  is  defined  as 

Cas  -   "J  *f  ■» 

or  for  constant  altitude  and  slow  speeds , 

AS  S   JV 

C^g  is  seen  to  be  proportional  to  the  negative  of  the  thrust  change 
with  respect  to  the  velocity  change.    Inequality  (3)  may  then  be  written  as 

Ineq.  (4)  1    - -    -77-       >0 

A  more  sophisticated  approach ,  based  on  the  suppression  of  altitude 
disturbances  by  the  pilot's  use  of  the  elevator ,  and  originating  with  the 
basic  equations  of  motion,  would  produce  essentially  the  same  inequality, 
with  a  small  modification  factor  included  „  (See  Refo  2)    The  resulting 
inequality  is 
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However ,  the  correction  term,  C     .  ,  is  of  relatively  small  magnitude 

D/a 

and  may  be  ignored.    From  this  same  analysis  the  time  constant,  "7"", 
of  the  subsidence  of  an  airspeed  error  is  shown  to  be  proportional  to: 

Typical  values  of  the  parameter  I  for  the  F-8  are  given  below. 
At  1.1  VgL/  I  =  .22  -  o548  =  -  .328,  showing  distinct  speed  instability. 
At  1 . 2  VgL,  I  =  .  200  -  .  222  =  -  .  022  ,  a  less  unstable  value . 

From  an  examination  of  the  speed  instability  parameter,  Ineq.  (4) , 
it  will  be  seen  that  positive  stability  can  be  produced  if  either  drag  in 
trimmed  conditions  (Cnr>)  or  C^o  is  artificially  increased. 


Ct>q  could  be  increased,  for  example,  by  extending  the  speed 
brakes  during  an  approach.    This  would  accomplish  two  objectives;  one, 
increase  C-qq  in  the  stability  formula,  and  two,  cause  the  engine  to 
operate  in  a  higher  RPM  range  to  produce  the  added  thrust  necessary, 
thus  effectively  causing  a  slight  decrease  in  the  inherent  engine  time 
constant. 

Remembering  the  definition  for  C»o  #  the  drag  coefficient  corre- 
sponding to  the  thrust/speed  variation,  it  will  be  seen  that  an  artificial 
means  of  increasing  thrust  with  decreasing  airspeed  accomplishes  the 
desired  objective .    This  is  exactly  what  automatic  throttle  or  approach 
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power  compensator  devices  are  intended  to  do.    Still  another  method  of 

increasing  C       is  by  the  use  of  an  anti-drag  technique.    This  may  be 
Ad 

accomplished  by  having  a  drag  producing  device  (such  as  side  extending 
speed  brakes)  extended  during  the  approach.    Then,  by  using  an  airspeed 
controller  type  device  to  regulate  the  amount  of  drag  produced  (by  opening 
or  closing  the  speed  brakes)  in  response  to  thrust  (or  anti-drag)  required, 
airspeed  deviations  can  be  closely  controlled „    This  particular  method, 
within  its  limitations  of  having  enough  anti-drag  available,  theoretically 
has  the  advantage  of  a  much  shorter  time  lag  than  an  engine  response 
scheme  o 


3.        Artificial  Speed  Stability  through  Automatic -Throttle  Techniques 

There  is  a  choice  of  inputs  to  an  automatic -throttle  type  controller, 
including  various  combinations  of  derivatives  and  integrals  of  airspeed 
error,  pitch  attitude  error,  angle  of  attack  error,  and  normal  acceleration. 
There  are  also  many  different  philosophies  by  manufacturers  concerning 
the  ideal  combination  for  best  performance  under  all  conditions  „    Two 
representative  systems ,  one  utilizing  angle  of  attack  and  normal  acceler- 
ation, and  the  other  using  airspeed  error  and  pitch  angle,  were  chosen  for 
investigation,, 

A)    Angle  of  Attack,  Normal  Acceleration  Inputs 
The  Specialties,  Inc.  Automatic  Power  Compensator,  (APC) ,  an 
auto-throttle  in  use  in  the  F-8  and  F-4  aircraft,  is  an  example  of  a  control 
using  these  inputs.    A  block  diagram,  typical  of  this  type  of  system  is 
shown  in  Fig.  3.    In  this  system,  constant  angle  of  attack  is  maintained 
in  the  landing  approach  using  the  two  sensed  inputs  of  instantaneous  angle 
of  attack  and  acceleration  normal  to  the  glide  slope  in  g  units  „    The  normal 
g  input  supplies  anticipation  whenever  the  aircraft  is  not  in  a  one  g  flight 
condition.    A  thrust  command  is  calculated  in  the  system  computer  depend- 
ent upon  error  signals  and  chosen  gain  and  temperature  factors.    This 
thrust  command  then  becomes  an  input  to  the  engine  which  in  turn  adds 
an  appropriate  thrust  to  the  airframe.    This  thrust,  acting  as  the  forcing 
function  to  the  airframe  modes  of  motion,  changes  the  velocity  of  the 
aircraft.    Thus  the  APC  computer  indirectly  controls  airframe  velocity 
through  the  equations 
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Eq.  (6)  ._     £V  =  ^4K   +K^nj 

where  the  K  values  are  appropriate  shaping,  filtering  and  simulation 

networks . 

Or,  looking  at  it  differently,  the  C       of  Ineq.   (4)  has  now  been 

Ao 

given  a  greater  value  f  Cfi&  3T  ~  s_T  J 

artifically,  thus  making  the  speed  stability  parameter  more  positive  or 
stable . 

The  transfer  function  equations  for  this  auto-throttle  as  supplied 
by  Specialties,  Inc. ,  are  also  shown  in  Fig.  3.    It  will  be  noted  that 
the  throttle  servo  is  simulated  by  a  .  1  second  time  constant  first  order 
lag.    Angle  of  attack  input  is  treated  by  both  an  integrator  and  lag  net- 
work.   Normal  g  input  is  also  integrated  „ 

The  results  of  the  analog  simulation  runs  using  this  type  of  control 
are  reported  in  Section  IV   of  this  paper,  Experimental  Methods  and 
Results  o 

B)    Pitch  Angle ,  Airspeed  Change  Input 

Flight  path  angle ,  7f  ,  is  given  by  the  equation 
Eq0   (7)  Y  =  0--<X 

During  constant  airspeed  flight,  pitch  attitude  change  (A#*)  and 
flight  path  angle  change  (Atf )  are  proportional;  but  dynamically,  A  &— 
occurs  before  AX .    Thus  the  means  of  anticipating  flight  path  angle 
changes  is  available  through  pitch  attitude  changes  recorded  by  a 


vertical  gyro  mounted  in  the  aircraft.    This  is  the  theory  underlying  the 
pitch  angle,  airspeed  change  type  of  automatic -throttle  controller,, 
Through  the  use  of  a  properly  designed  network,  these  inputs  make  it 
possible  to  minimize  transient  airspeed  changes  and  to  eliminate  steady 
state  airspeed  changes  resulting  from  variations  in  flight  path. 

The  governing  equation  for  this  controller  may  be  of  the  form: 

Here  the  K      term  is  inserted  as  a  steady  state  error  washout  term. 
The  lag  term  (  /  +  Ja^-)  simulates  the  airspeed  sensor  and  filter  „    The 
filter  attenuates  the  short  period  wind  gust  signals  entering  the  system 
and  acts  as  a  smoothing  influence. 

This  equation  was  reduced  to  a  linear  open  loop  transfer  function 
of  the  auto-throttle  controller.    The  block  diagram  for  this  system s  com- 
bined with  the  analytically  derived  block  diagram  and  transfer  functions 
of  the  airframe  and  engine,  are  shown  in  Fig.  40    The  stability  derivatives 
for  the  F-8  in  landing  configuration  were  used  for  the  airframe  transfer 
functions .    The  derivation  is  shown  in  Appendix  I . 

The  system  open  loop  transfer  function  was  calculated  as: 


Eq.  (9) 


**      nY-^  (£s\J£-)j£\r 
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where;      l^j/\r  "  J[(,  +/./^)(,  *-.  /<) 


fel 


uA     _  .£,2**10-*  fi  +  3724.  +  !  8.  ?4>  +  £3.154} ) 
/AT)         -       J^ 


This  equation  was  programmed  on  the  digital  computer  to  produce  a 
series  of  Nyquist  plots ,  each  utilizing  different  gain  constants  . 
Fig.  5  shows  the  Nyquist  plot  for  the  equation  using  the  same  gain 
constants  as  were  used  in  the  analog  simulation .    From  this  figure  the 
open  loop  gain  margin  is  seen  to  be   19. 2 . 

The  system  closed  loop  transfer  function,  —j-rr-     )  was  then  derived 
This  is  shown  symbolically  as  follows: 


„    ,in,       u.     _    IallJat(atJaf 

Eq.   (10)      --—       = 


[allJat{&T; 


^*~   >i^lr(&m^<U^. 


A  digital  program  was  devised  to  calculate  the  values  of  the  poles 
of  the  closed  loop  transfer  function  for  various  values  of  auto -throttle  K's 
This  procedure  would  have  produced  families  of  Root  Loci  Plots  „    How- 
ever, because  of  lack  of  time  and  computer  malfunctions,  this  particular 
investigation  was  not  completed.    Equation  (10)  would  be  used  in  an 
optimizing  analysis  for  gain  constants  „ 

In  addition  the  equations  of  the  auto-throttle  were  set  up  on  the 
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the  analog  computer.    This  was  incorporated  into  the  complete  system 
using  approximately  optimum  values  of  K^.     0  Ku,  and  Kyu»    All  simu- 
lation was  done  in  real  time.    These  results  are  also  reported  in  Sec- 
tion IV  .     Analog  mechanization  is  shown  in  Appendix  III  .    A  composite 
block  diagram  showing  the  relative  positions  of  the  two  automatic 
throttle  systems  in  the  airspeed  loop  is  presented  as  Fig.  6. 
C)  Auto-Throttle  Controller  Linked  to  Anti-Thrust  Device 
Another  possible  approach  to  the  airspeed  instability  problem  is 
the  use  of  an  auto-throttle  controlling  an  anti-thrust  or  drag  device. 
This  could  be  accomplished  byflyingthe  approach  with  speed  brakes 
extended.    The  output  of  the  auto -throttle  would  then  reduce  or  in- 
crease the  drag  by  closing  or  opening  the  speed  brakes  as  required. 
Thus ,  instead  of  supplying  additional  thrust  when  the  airspeed  drops , 
the  speed  brakes  would  be  closed  a  proportionate  amount,  reducing  the 
drag  and  in  effect  accomplishing  the  same  purpose  as  an  increase  in 
thrust.    Note  that  the  time  lag  of  the  hydraulically  actuated  speed  brake 
system  is  much  less  than  that  of  the  engine,  thus  allowing  a  taster  res- 
ponse.    The  time  constant  of  the  hydraulic  system  is  of  the  order  of  „  1 
second,  while  that  of  the  F-8  engine  is  about  1.16  seconds  „    This 
approach  was  also  investigated  on  the  analog  computer  by  replacing  the 
engine  with  a  first  order  time  lag  of  .  1  seconds  (simulating  the  hydraulic 
system)  and  using  the  thrust  command  output,  of  the  auto -throttle  com- 
puter to  control  the  speed  brakes . 


Obviously  the  limiting  feature  to  such  a  system  is  the  amount  of 
drag  the  speed  brakes  are  capable  of  supplying  „    Of  course     it  goes 
without  saying     that  such  a  device  can  not  be  installed  on  an  aircraft 
incapable  of  completing  an  approach  while  varying  the  positioning  of 
the  speed  brakes  „    An  airplane  in  this  category  is  the  F-8  which  has  a 
lower  fuselage  speed  brake.    Insufficient  deck  clearance  is  available 
when  this  brake  is  extended ,  thereby  prohibiting  its  use  during  an 
approach . 
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4.         Experimental  Methods  and  Results 

A)    Equations  and  Analysis 

The  basic  3-degrees  of  freedom  aircraft  longitudinal  equations  of 
motion  (in  British  notation)  used  in  the  analysis  are  shown  below .    The 
equations  are  based  on  the  small  perturbation  linear  theory.    Lower  case 
letters  represent  perturbation  quantities.    Upper  case,  subscripted  o, 

letters  represent  initial  steady  state  conditions  .    Body  axes  are  used 

^ r  . 

throughout.    The  assumptions  made  are  as  follows: 

1.  The  airplane  has  a  longitudinal  plane  of  symmetry. 

2.  The  direction  of  the  relative  wind,  in  steady  flight,  lies  in 
the  plane  of  symmetry,  and  in  steady  state,  all  angular  velocities  are 
zero. 

3.  Initial  flight  condition  is  wings  level. 

4.  Products  and  squares  of  perturbation  velocities  are  small  and 
are  ignored,  and  sines  of  all  perturbation  angles  are  approximated  by 
the  angle  itself  in  radians. 

5.  The  longitudinal  modes  of  motion  are  independent  of  the 
lateral  modes . 

6 .  Structural  deformations  are  not  considered . 

Three  degree  of  freedom  non-dimensional  equations  of  motion, 
British  notation:  (See  Ref .  3)  Figure  7  displays  angular  definitions  and 
signs. 
X  Force  Equation: 

Lt  ^-]£  -*-  g  - O  -£-&]» = K&T 
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Z  Force  Equation: 
Moment  Equation: 

L   4  .lie  L  a  ^7-  -cr^h/o  LW-  ZotrJ   1  <•<>  1  l  i<>  J 

where  symbols  are  defined  as  follows: 
"t      =  07  re    »    ^^  sees  . 

0L7-        dt 

*>*    =  jt_       real  time 
J  £        air  sees 

u       =  perturbation  horizontal  velocity,  ft./sec<, 

w      =  perturbation  vertical  velocity ,  ft. /sec. 

*}      =  elevator  deflection,  positive  down,  radians 

A  /     =  change  in  thrust  required  ,  lbs  0 

OC     „  angle  of  attack,  radians 


<6L 

Vo 


$■    =  pitch  angle  ,  radians 

Juu 

CZ~a~i  v  Inertia  Coefficient 


U    = 


T  z 

Lo  =  Moment  of  inertia,  slug  -  ft 

Inn    =  mass  ,  slugs 

Jtt  =  tail  length,  ft. 

/*\  =  pc*A        =  relative  density 

C  =  mean  wing  chord,  ft„ 

S  =  reference  Wing  area,  square  ft. 

j(_  =  Force  parallel  to  x  -  axis,  lbs. 

£    =  Force  parallel  to  2  —  axis  ,  lbs  . 

j*|_  =  Moment  about  lateral  axis,  lbs.  -  ft, 
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y-u     = 


mur~ 


distance  from  thrust  line  to  CG,  ft„ 

1/2     {-tflfcj       non-dimensional  Force  Stability  derivative 


(Ll\  an 

(-U7  x, 


and  similar  identities  obtained  by  permuting 
2    and  \i„ur  independently  o 


x0     =  1/2     l^97f?^/     non-dimensional  stability  derivative  and 

similar  identities  obtained  by  permuting  x, 

2- 

xT     =  1/2    A  pn-T? )      non-dimensional  stability  derivative  and 

JfLI similar  identify  obtained  by  permuting  x 

^ T  and  ^_  „  in  and   T 

#\ 

25cJ  non-dimensional  moment  stability  derivative 

£t  \fj^\  and  similar  identities  obtained  by  permuting 

^VV  u  and    U-r 

x~.    =  l/2&(r£flj2s}       non-dimensional  stability  derivative  and 

s  &-  similar  identities  obtained  by  permuting  x, 

and    /X^- 


^w^ 


1/2=.  UyJ'StJ     non-dimensional  stability  derivative 

mT    =  y  -Jl  non-dimensional  stability  derivative*  h  = 

'  j&  distance  from  thrust  line  to  CoG=  position 

The  necessary  stability  derivatives  were  evaluated  by  making  use 
of  a  longitudinal  dynamic  stability  Fortran  program,  designated  LONGSTAB„ 
The  program  was  compiled  while  undertaking  the  course  in  dynamic  air- 
craft stability  given  at  the  Postgraduate  School .    This  program  computed 
the  necessary  British  stability  derivatives  and  then  used  these  derivatives 
in  the  aerodynamic  equations  of  motion.    The  resulting  stability  quartic 
was  then  solved ,  by  the  program ,  for  the  phugoid  and  short  period 
complex  roots.    Finally,  the  periods  and  times  to  damp  to  1/2  amplitude 


were  computed  for  both  modes.    Necessary  aircraft  parameters  used  by 

the  program  were  taken  from  F-8  data  supplied  by  Ling-Temco-Vought ,  Inc. 

or  calculated.    The  program -computed  stability-derivatives  used  in  the 

investigation  are  shown  in  Table  I,  as  are  the  other  aircraft  parameters. 

All  program-computed  derivatives,  after  necessary  conversions  were 

found  to  be  in  close  agreement  with  those  listed  in  Ref .  6    „    The  computer 

program  is  explained  in  Appendix  II.    A  sample  printout  is  shown  in 

Table  II. 

From  the  program  printout,  the  phugoid  period  is  38 . 5  seconds  and 

that  of  the  short  period  mode  is  6.  OS  seconds.    Times  to  damp  to  1/2 

amplitude  are  29.9  and  1.73  seconds  respectively.    Undamped  natural 

frequency.  U)      L        .j«  is  calculated  as  . 189  radians/sec. ,  X  as  .123. 
^n  phugoid  ° 

Short  period  Co    is  calculated  as  1  „  19  rad./sec.  ,  T  as  .360. 
n  J 

Ref.  6  lists  a  phugoid  (J     of  0.188  rad.,/sec.  ,  and  a   "Tof  .12. 

n  ^ 

Short  period  results  found  were  listed  as£j    =1.19  rad./sec. ,  ^  = 

n 

0.35.    The  results  given  by  LONGSTAB  and  those  given  by  Ref.  6  are  in 
excellent  agreement. 

Thus,  as  previously  remarked,  it  is  seen  that  the  phugoid  is  very 
lightly  damped. 


Ref  0  6  reports  on  the  results  of  an  investigation  conducted  along 
similar  lines  as  the  subject  of  this  paper.    The  two  investigations  were 
performed  concurrently  and  independently  of  each  other.    After  com- 
pletion, the  results  of  the  author's  analysis  were  verified  insofar  as 
possible  with  those  given  in  Ref.  6. 
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These  non-dimensional  equations  of  motion  were  dimensional! zed 
for  the  analog  simulation  ,    The  engine  time  lag  and  auto-throttle  blocks 
were  next  added  „    Real  time  stick  and  throttle  simulators  were  then  in- 
troduced into  the  system ,  together  with  the  F-8  linearized  longitudinal 
control  system  dynamics  supplied  by  Ling-Temco-Vought.    All  analog 
runs  were  conducted  in  real  time.    The  analog  computers  used  were 
Donner  3100  and  3400'  machines.    An  eight  channel  Brush  recorder  was 
used  to  record  the  time  histories  , 

The  following  types  of  computer  runs  were  then  conducted: 

(a)  Transient  analysis  runs  using  step  inputs  of  elevator  deflection, 
horizontal  velocity  gusts ,  and  pitch  attitude;  first  with  basic  airframe 
alone  and  then  with  the  two  basic  types  of  auto-throttle  controllers  in- 
corporated o    Finally  sample  runs  were  made  with  the  speed  brake  device 
replacing  the  engine  „ 

(b)  Optimizing  runs  using  triangular  inputs  of  elevator  deflection, 
to  determine  the  optimum  gain  for  this  type  of  forcing  function, 

(c)  Sinusoidal  gust  inputs, 

(d)  Real  time  pilot,  controlled  runs  using  a  simulated  glide  slope, 
with  and  without  the  aid  of  the  auto-throttle.    Horizontal  turbulence  was 
also  added  to  determine  the  controller  effectiveness. 

Analog  mechanization  diagrams  are  shown  in  Appendix  III ,    Poten- 
tiometer settings  are  shown  in  Table  1  to  Appendix  III.    Time  histories  for 
these  runs  are  displayed  in  Figs.  10  to  35. 
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B)    Evaluations  of  Experimental  Analog  Runs 

The  systems  referred  to  are  as  follows; 

System  1 ,  Inputs  of  angle  of  attack,  and  normal  acceleration 
in  g  units  . 

System  Equation:       AT=-L_  f  .****    -  KavAV      Kx  A*l 

System  1-a  gain  constants? 
K^     =     1970  lbs. /degree 
K^v  =     280    lbs ./knot 
K^     =     183    lbs  ./degree  -  sec. 

System  1-b  gain  constants; 
K^    =     3420  lbs  ./degree 
KAV   =     72      lbs. /knot 
K-^     =     190    lbs. /degree  -  sec. 

System  2 ,  Inputs  of  pitch  angle  „  airspeed  error 

System  Equation:     AT  -  K&A&- -  -~J\    '  ^  +  ^H±l 

Gain  constants; 

Kp.    =     318  lbs. /degree 

Kn    =     360  lbs. /knot 

Kjn,   =     22.3  lbs. /knot  -  sec. 

1)    Step  input  of  1.32°  up  elevator  deflection 

Fig.  10  is  the  analog  time  history  response  to  a  1.3     step  of  up 

elevator  deflection ,  using  the  System  l~a  auto -throttle.    Airspeed  dis= 

plays  a  slight  oscillatory  tendency,  reaching  an  equilibrium  speed  of 

+  .6%  of  V     (.8  kts .  above  V_) .    Thrust  is  applied  smoothly  in  response 
o  ° 

to  demand. 

Fig.  11  shows  the  response  of  System  l-b„    This  system  had  the 
gain  factors  optimized  for  a  step  input  of  horizontal  gust  velocity.    Note 
the  high  airspeed  overshoot  due  to  step  elevator  input. 
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Fig.  12  is  the  same  situation  usii      S    s         ?,.    Airspeed  is  muc 
more  highly  damped  as  compared  to  F  3      Equilibrium  airspeed  is 

-l%ofV     (1.39  lets,  below  V0)o    Thrust  is    .     o  applied  smoothly. 

2)  5  knots  of  horizontal  gu  31 

Fig.  13  displays  the  record  of  i  response  to  this 

input.    The  system  is  oscillatory  v.        n    appro*  n  ate  'lamping  ratio  of 
$  ■  .15,6^   -  .  677  rad/sec.    The  e  oc        o  ershoot  is  61%  of 

initial  error.    These  values  of  g  a  e   Less  than  optimum  for 

this  type  of  input  disturbance . 

Fig.  14  displays  the  response  of  .    stem         r  which  has  teen 
optimized  to  this  input  disturbance  c    I         alent  t   is  now  „58,  £j    is 
.73  rad/sec.    Thrust  response  is  somewhat  oscillatory .    Final  value  of 
airspeed  is  +.  2%  of  V   . 

Fig.  15  is  the  record  oi:  Syste."  r  ows  imped  response 

with  final  airspeed  of  about  -2     °A  .  changes  are  less  than 

that  required  by  System  1 . 

3)  4     step  input  of  pitxn     /     Le0 

Fig.  16  shows  the  response  o*  ,       :e  to  this  input.    Initial 

velocity  transient  is  +.8%  V    or  .    S"  e  tate  error  is  +.4%  V    . 

o  o 

Thrust  application  is  slightly  osci 

Fig.  17  shows  System  1-fc  respoi  se.     T-  itial  velocity  transient  is 

.6%V    .    Steady  state  error  is  also 
o 

Fig.  18  displays  the  response  o    S    stem  -       [nitial  velocity  transient 
amounts  to  +1.6%  V    .    Airspeed      m  >  s  to  +  A%  VQ  steady  state  error  in  19 


seconds,  reaching  less  than  1%  error  in  9  seconds. 

Fig.   19  shows  the  response  of  the  basic  airframe  to  this  same 

input  disturbance,  without  the  aid  of  the  auto-throttle.    Airspeed 

oscillation  is  +1%  V    .  with  a  steady  state  value  of  about  ,  i%  V    , 

o  o 

4)    Optimizing  runs  using  triangular  elevator  input 

Fig.  20  shows  System  1-a  under  the  influence  of  a  .02  ops 

triangular  elevator  deflection  (maximum  amplitude  I  ,  32  )  forcing 

function.    Airspeed  error  is  kept  within  a  range  of  +2%  to  •=•!%  of  V    . 

Thrust  varies  within  the  range  of  +1600  to  =2400  lbs*  about  the  trim 

setting.    Thrust  application  follows  elevator  movement  closely. 

Fig.  21  shows  the  reaction  of  System  1-b.    Note  how  airspeed 

fluctuates  +  4%  V     (+5,56  kts.) .    Thrust  has  very  noticeable  oscillatory 
o 


tendencies  and  gyrates  within  rather  wide 

Fig.  22  shows  System  2  under  the  influence  of  this  same  forcing 
function .    Airspeed  is  seen  to  vary  +  i  ,  7  kts  ,  or  +  1 = 2%  of  VQ .    Thrust 
varies  approximately  2400  lbs,  about  trim  thrust.    Application  of  thrust, 
however,  is  smoother  and  displays  less  hunting  characteristics  than 
shown  by  System  1-a. 

Figs.  23 ,  24  and  25  show  the  response  of  these  same  systems  to 
a  triangular  elevator  forcing  function  of  the  same  magnitude  as  previously 
but  of  .06  cps  frequency. 

Fig.  26  shows  the  reaction  of  the  airplane  to  this  same  fore  in  g 
function  without  the  benefit  of  the  auto -throttle  „    Note  airspeed  deviation 
is  less  than  that  obtained  with  System  1-b  engaged. 


5)  Human  pilot  controlled  flight  on  simulated  glide  slope 

For  this  series  of  runs  a  pilot  flew  the  simulated  airplane  on  the 
glide  slope,  using  the  auto-throttle <>    Glide  slope  reference  was  obtained 

from  a  zero  centered  panel  meter.    The  zero  voltage  level  was  considered 
the  correct  glide  slope.    This  is  shown  on  Channel  1  of  Figs.  27,  28  and 
29 .    The  pilot  was  then  told  to  manipulate  the  stick  to  fly  the  aircraft  to 
a  new  4    higher  glide  slope*    This  corresponded  to  a  +•  3  4  volt  reading  on 
channel  1 ,  or  14  volts  on  his  reference  meter.    This  simulated  a  low 
"meatball"  on  a  mirror  approach  and  the  flight  path  correction  necessary. 
This  new  flight  path  angle  was  held  momentarily.    The  pilot  was  then  told 
to  fly  the  airplane  down  to  a  -  4     glide  slope  ,  after  which  he  was  to 
return  to  the  normal  flight  path  (0  volts) . 

The  response  of  the  aircraft,  and  speed  control  is  shown  in  Figs. 
27,  28  and  29  for  Systems  1-a,  1-b  and  2  respectively. 

It  will  be  noted  that  Systems  1-a  -md  2  maintain  airspeed  close  to 
within  1%  of  VQ/  whereas  System  1-b  has  much  greater  airspeed  flue: 
tions.    Thrust  application  is  again  much  more  oscillatory  in  System  1-b. 

6)  Human  pilot  controlled  flj  jht  on  glide  slope  with  sinusoidal 
input  of  5  knots  horizontal  yelo< 

Figs.  30  and  31  are  records  of  System  1  runs  ul  .  ...if     a  pilot 
controlled  stick  with  a  forcing  function  sinusoidal  input  of  5  knots  of 
horizontal  velocity,  frequency  .05  cps.    The  pilot  wa        .  tin  attempting 
to  fly  the  glide  slope  4    up  and  down. 

Fig.  32  is  the  record  of  the  response  to  this  same  forcing  function 
under  the  influence  of  the  pilot  and  the  System,  2  auto-throttle. 
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Fig.  33  displays  the  record  of  the  pilot,  attempting  to  maintain 
original  glide  slope  while  the  aircraft  was  being  influenced  by  a  sinu- 
soidal u  turbulence  forcing  function,  without  the  benefit  of  the  auto- 
throttle  . 

7)    Speed  brake  anti-thrust  system 

Fig.  34  shows  the  history  of  the  speed  brake  anti -thrust  system 
in  response  to  a  5  knot  horizontal  velocity  step  input „    Mechanically, 
the  System  2  auto-throttle  was  used  to  supply  a  thrust  command,  but 
this  thrust  command  was  used  to  drive  a  ,1  second  time  constant  anti- 
drag  device,  simulating  the  movement  of  the  speed  brakes.    It  will  be 
noted  that  in  comparison  with  Fig.  15,  (standard  System  2  setup)  response 
and  damping  is  faster.    It  required  9  seconds  for  the  transient  to  steady 
within  1%  of  V    in  the  case  of  Fig.  34,  while  it  required  11  seconds  for 
System  2. 

Fig.  35  shows  the  response  of  the  anti -thrust  system  to  the  same 
forcing  function,  but  utilizing  different  gain  factors,  Response  is  more 
oscillatory . 

C)    Summary  of  Analog  Results 

From  the  foregoing  evaluations  it  can  be  seen  that  the  response  of 
the  different  systems  investigated  is  dependent  upon  the  sensed  variables 
and  the  applied  forcing  functions.    System  l~a,  optimized  to  an  elevator 
deflection  input,  is  marginal  to  unsatisfactory  in  performance  when 
subjected  to  a  horizontal  gust  input.    On  the  other  hand,  system  1-b, 
using  the  same  sensed  variables,  but  with  gains  optimized  to  horizontal 
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velocity  input,  does  not  perform  well  when  subjected  to  elevator 
deflection  inputs  .    System  2  ,  utilizing  different  sensed  inputs  „  also 
displays  varying  reactions  to  each  of  the  forcing  functions . 

In  general,  it  can  be  seen  that  System  2  has  the  capability  of  a 
smoother  application  of  thrust  in  response  to  an  airspeed  error „    System  '. 
has  the  capability  of  a  more  highly  damped  response,  but  causes  a  more 
erratic  application  of  thrust.    Further  investigation  as  to  the  optimum 
values  of  gain  constants  is  required  to  accomplish  the  desired  task  more 
efficiently , 

It  is  further  recommended  that  additional  sensed  variables  or 
integrals  of  such  terms  be  added  to  the  analog  set-up  to  help  optimize 
system  response.    A  fruitful  area  might  be  to  incorporate  a  normal  g 
input  to  System  2  to  help  the  damping  characteristics  „    It  is  anticipated 
that  further  investigations  along  this  line  will  be  conducted  at  the  Uo  S„ 
Naval  Postgraduate  School.    These  investigations  will  attempt  to  utilize 
the  simplified  equations  obtained  by  the  elimination  of  the  short  period 
mode,  to  develop  a  more  tractable  analysis 0 
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5.         Conclusions 

The  two  major  approach  power  compensator  systems  investigated 
appear  to  have  the  capability  of  performing  the  desired  task.    The  major- 
difference  are  degree  of  anticipation  and  erratic  ism  of  thrust  application. 
With  proper  choice  of  gains  and  input  variables ,  it  seems  feasible  that  a 
composite  system  can  be  designed  to  provide  the  advantages  of  both 
individual  systems. 

The  speed  brake  or  anti-thrust  system  has  the  capability  of  pro- 
viding faster  response,  but  is  limited  by  drag  area  available  and  the 
physical  feasibility  of  approaching  with  speed  brakes  extended. 

The  preceding  investigation  has  firmly  convinced  the  author  that 
an  auto -throttle  mechanism  is  an  important  contribution  to  carrier  aviation 
safety.    An  aircraft  that  is  difficult  to  fly  in  the  approach  because  of 
inherent  speed  instability  problems,  c  u    r    ve  its  whole  character  changed 
by  an  auto-throttle .    By  removing  the  undeslred  instability  an  auto-throttle 
relieves  the  pilot  of  concentrating  on  airspeed  control,  thus  allowing  more 
precision  in  glide  slope  and  line-up  control.    This  advantage  by  itself, 
should  tend  to  reduce  the  frequency  of  ramp  accidents  „ 

In  addition,  when  an  automatic  landing  system  is  considered,  the 
pilot-controlled  airspeed  loop  must  be  replaced  by  an  automatic  system 
to  provide  consistency  and  predictability  of  response. 

Still  another  advantage  of  the  automatic  airspeed  control,  system  is 
the  reduction  of  minimum  approach  speed  below  that  recommended  by  the 
pilots  when  airspeed  is  controlled  by  manual  throttle  manipulation.    In 
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addition  the  automatic  device  can  be  de      jne     to  pro  ride  f  cipation 

and  instant  response  that  can  only    >e  obtained  from  actual  pilot 

experience.    Thus  it  is  seen  that  the  auto         ottle  <  ct  as  the  great 

equalizer  in  carrier  approach  trainii 
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APPENDIX  I 


Derivation  of  Transfer  Functions 

The  dimensional  derivatives  shown  in  Table  I  were  used  'o  calculate 
the  analytic  transfer  functions.    Calculated  values  of  these  derivatives 
using  F8  data ,  are  shown  in  Table  II. 

The  longitudinal  Equations  of  Motion  written  in  Matrix  Laplace  Form 
are  (neglecting  negligible  terms); 

<i-X*    -Xv,    +?  \/'Ltd4,\     (°    X*T\    ft*\ 


Mn.  -(iMv+M^  4-l\y\&C4)j 


The  Transfer  function  /~t~1  Airframe  then  can  be  written  as 
(in  determinant  form); 

0  "Xur  M 

2je       4-^Xw  -4U0 

Hie-    -(aMw-^Mw)     4-l^fA. 


hi 


Lc 


D, 


Di 


Where  D    =  determinant  of  homogeneous  equations 


Expanded  D    =  A-4.4  +Ba    +  C*?  ^L^  +  E 
Where  A  =  1 

B  =  -  (tAf  +Xlc  +?ur  +U>Mur) 


Expanded       Mu-  ~    Dec  A.     +   CU  -d.    +  D<^- 
Where      Bu  =     £jeXocr 

Substituting  values  for  these  derivatives 


N 


u. 


By  a  similar  procedure  , 

(&\  _    -^.SqC^1  +-.2/554  +  .OQS68] 

V<£eMiHVi*»e.  D« 

(Positive  elevator  deflection  is  down) 

The  determinant  for  calculating  1—^] 


can    e  written  as 


Expanded: 


(& 


Wirfteme 


D, 


Where 


Ar=       X 


AT 


BT  =      -X*t[  Zur  +M  J  J  +  T/o  M(Xr 

CT=       XAT^Mj-UoM^I-MATlUXur-g] 

DT  =       AUt  3  ?ur 

After  substituting  in  values:  v       ',l  ^ 

>Qol*5  £a  f  .g043«-  +I/S14.  T-.0430J 
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APPENDIX  -X 

TABLE     J-A 
LONGITUDINAL  DIMENSIONAL  STABILITY  DERIVATIVE 
PARAMETERS 
(STABILITY  AXIS  SYSTEM) 


QUANTITY 

IN  TERMS  OF  BASIC  STABILITY  DERIVATIVES 

IN  TERMS  OF  NON- 
DIMENSrONAL 
STABILITY  DERIVATIVE 
PARAMETERS 

DIMENSIONAL 

NON-DIMENSIONAL 

definitions 

UNIT 

X. 

J  1  ;*x 

1 
sec 

■«f!  (-c0-cDb) 

Ix. 

T 

x. 

1 
sec 

■S^-Si) 

1 

x«. 

.1  M. 

■  3SB 

ft 
scc2rad 

■f2(-c".) 

r    '■ 

*        2. 

.12? 

■  3u 

j        sec 

•f^-O 

r 

\ 

•IS 

i 

.     _L 

^        sec 

'^(-<VCo) 

M. 

z* 

.1  23 

■  air 

1 
1 

"«(S) 

■** 

2, 

a   cq 

ft 

•^(S) 

'%** 

sec-rad 

Z<« 

.1  3S. 

m  as. 

1 

2"       ^        '  «£,) 

.uZl 

sec'rad 

.JL2a 

I,   eu 

1 

•*if(c-*c-.) 

•  -«.  a 

TC       ■ 

sec- ft 

M. 

,_L  » 

1 

sec- ft 

■^c-« 

TC       " 

«. 

1    W 

"  I,  s« 

1 
ft 

■^s* 

*rDn% 

M, 

I,   2<I 

1 

41r      '- 

-I- 

T     * 

sec  rad 

"*r 

JL    r"' 

1 

.  *f£s  c 

re     *g 

sec2rad 

l> 
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TABLE  J-& 

LONGITUDINAL  NON-DIMENSIONAL  STABILITY  DERIVATIVES 
(STABILITY  AXIS  SYSTEM) 


BASIC  NON-Pr^ENSIONAL  STABILITY  DERIVATIVES 

NON-DIM  ENSiONAL 

STABILITY  DERIVATIVE 

PARAMETERS 

TOTAL  A IH FRAME 

THEORETICAL  HORIZONTAL       . 
TAIL  CONTRIBUTION 

DEFINITIONS 

UNIT 

r        DRAO 

1 
1 

I 
1 

».'  (Cd-C^) 

ra3 

"-•|<cu-CDa) 

C        .-5 

_1 

rad 

Ik     2      o,e 

fc        OS 

1 

c     -y  ** 
••a    ?    atl 

_L 
1 

%  ■  (<V<:Lj 

1 
rad 

N."O*H0 

*---z-<-Cta-cn) 

Li  »S) 

I 

rad 

rc  i  .  oc     1»  *S  h.  *i 
LL<iJ„  "  l<rh  q  s  '  M 

tL.'J-    C.    . 

•      4       ta 

C    .*£i. 

1 

rad 

LSL  "Si «  «  ■ 

'«?% 

1 
rad 

.       L',     u«h  a   s  as; 

-i- iS. 

c.  •  -Jj- 

qSe 

1 
1 

C      .   S    -5 
"a       2    »U 

1 
1 

■••ilvj'^.-V 

e    •** 

'a     a« 

1 

i  ad 

N."=N. 

■••i(^)'c-. 

""""•si 
"MS) 

1 

rad 

1 

r<td 

"••t(^)'s 

■«r.   >•• 

1 

rad 

c. E  c. 

V*(«N 
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Appefldiy  X    TzbleJL 
*       4-X^  4    ^  74**0}"  5 


2*.    ^flTSf-a-O  =  (£og)  (-2SW-H7)-      -.+16    -J- 

Xc,      -?WS(-Cd^-QO  =  2°8  (-.117)        ~      -,4Ta  _L 
~TPT  ~ggj  sec 

Xur-  r  f^5  (C-Cdk)   =   ^^8    (.*70-Ui)-  -  .0I4ISJL 

X4T~  Cos  (j*&)         t=  DOITS'       4± 

Mat   -    -J±  r  -4.3-5-  *I0~6_L. 
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APPENDIX  II 
DYNAMIC  STABILITY  FORTRAN  PROGRAM 

Longitudinal  -  Program  LONGSTAB 

This  program,  once  entered  with  required  airplane  parameters , 
computes  the  longitudinal  stick  fixed  and/or  stick  free  (when  pertinent) 
stability  derivatives ,  and  the  corresponding  stability  equation  coeffi- 
cients.   It  then  solves  these  equations  for  roots,  periods,  and  times  to 
damp  to  1/2  amplitude  for  both  the  phugoid  and  short  period  modes .    In 
addition  the  value  of  Routh's  discriminant  is  determined.    The  sequence 
of  computations  can  be  repeated  for  as  many  different  flight  regimes  as 
desired. 

There  are  a  total  of  five  data  cards  required  for  the  stick-fixed 
only  solutions ,  with  an  additional  two  cards  required  for  a  stick-free 
solution.    Thus  for  a  combination  stick-fixed  and  free  solution,  there 
will  be  a  total  of  seven  data  cards  required  for  each  flight  regime 
desired. 

The  initial  data  card  controls  the  types  of  solutions  desired 
(stick-fixed,  free,  or  both),  and  the  number  of  consecutive  runs  to  be 
made.    A  list  of  program  symbols  with  their  meanings  appear  in  Table  I 
to  this  Appendix . 

The  second  card  contains  an  alpha-numeric  run  identification  or 
an  arbitrary  title . 

The  third  card  is  the  first  of  the  general  input  data  cards .    This 
card  should  have  entered  on  it: 
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St(tail  area),  Sowing  area),  CL     t,  "^taii*  SB  (body  area),  1 
(body  length) ,  f7  ,   g (acceleration) . 

The  parameters  are  entered  eight  to  a  card,  each  being  allowed  10 
columns.    There  are  80  columns  on  each  card»    Each  parameter  may 
be  entered  anywhere  in  its  assigned  ten  columns,  but  must  have  a 
decimal  point  included  in  it.    This  same  format  is  followed  for  each 
of  the  succeeding  cards  „ 

The  fourth  card  should  have  the  following  aircraft  parameters 
entered  on  it: 

W  (aircraft  weight,  lbs„);  —=•    ;   0o    in  degrees;  h,  (distance  from 

2 
e.g.  to  thrust  line,  ft„);  C    (coefficient  of  thrust);  C        ,°  Iw  (slug  ft  .); 

t  m(X       1X 

V  (forward  velocity  ft. /sec.) . 

The  fifth  card  should  have  the  following: 
CD;  CLc^;  eJTA  ;  CL  wing;  c  (wing  chord);  lt  (e.g.  to  a.c.  tail,  ft.); 
Xc.g.  (x  distance  of  e.g.  ,  ft.);  XAC  (x  distance  of  wing  a.c.  in  ft.) . 

Cards  2  through  5  are  repeated  for  each  different  run ,  stick  fixed . 
If  it  is  desired  to  compute  a  stick  free  solution  simultaneously,  then 
two  additional  data  cards  (6  and  7)  are  required  for  each  run. 

Card  6  contains  the  following  parameters: 

2 
CT  (tail  chord,  ft.);  K      (radius  of  gyration,  squared);^*-;  B   ;  B   ;    B' ; 
*  e  12        2 

z^ ;  z^    . 

Card  7  contains  the  following  parameters: 

2 
nu      ;  nu    ;  m     (mass  elevator,  slugs);  S     (area  elevator,  ft.). 
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Thus  for  one  aircraft  flight  regime  (say  Sea  Level ,  Mach  0.6), 
for  both  stick  fixed  and  free,  there  would  be  required  data  cards 
1  through  7.    For  each  succeeding  run,  cards  2  through  7  would  have 
to  be  repeated  with  the  new  parameters  entered.    Card  one  would  have 
to  be  suitably  prepared  to  reflect  the  total  number  of  runs  and  mode  of 
runs  desired.    See  Fig     1  this  Appendix  for  sample  data  cards  „    Note 
the  first  card  is  of  format  13  for  the  number  of  runs ,  and  II  for  mode  of 
operation.    All  other  data  cards  are  of  8F10.0  format  (each  card  holds 
8  fields  of  10  columns,  each,  one  data  number  to  be  entered  per  field, 
in  floating  point  format;  i.e.  ,  with  a  decimal  point  somewhere  in  the 
number.) 

The  program  prints  out  the  calculated  stability  derivatives  as 
well  as  the  computed  solutions  to  equations . 
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TABLE  I 


Program 
mnenonic 


PROGRAM  SYMBOLS  AND  MEANINGS 


Program  LONGSTAB 


Meaning 


NOS 

MODE 

ST 

SW 

CLALFAT 

ETATAIL 

SB 

BL 

RHO 

G 

W 

EPSALFA 

THETA1 

YH 

CT 

CMALFA 

YI 

V 

CSUBD 


Number  of  runs  to  be  made.    13  format,  card  1 
Decision  stick  fixed  (1),  or  stick  free  (2).    II.  card  1 


St,  tail  area,  card  3. 


S     ,  wing  area,    card  3. 


CL  .       ,  Tail  lift  curve  slope  =    card  3. 

^T    ,  Tail  efficiency,  CL/CL.    card  3. 
Sg ,  Body  area       card  3 . 
1R,  Body  length        card  3„ 

P    ,  Atmospheric  density,        card  3. 
g,  Acceleration  due  to  gravity,      card  3. 
w,  aircraft  weight,    card  4. 
^Ss      card  4. 

Q-Q,    ,  initial  angle  of  pitch,  degrees,    card  4. 
h,  distance  from  e.g.  to  thrust  line       card  4„ 


C.  ,  coefficient  of  thrust,      card  4 


mo^ 


per  radian .      card  4 . 


I      ,  moment  of  inertia .      card  4 . 

V,  forward  velocity,  ft.  /sec  „      card  4 


Cp.,  coefficient  of  drag=      card  5 
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CLALFA  Ct      ,  winglift  curve  slope,      card  5. 

EPIA  £TfA,      card  5 . 

CSUBL  CLwing,      card  5. 

C  ce  wing  chord,      card  5. 

TL  1. ,  e.g.  to  a.c.  tail,      card  5<> 

XCG  x  ,  x  distance  of  e.g.      card  5. 

XAC  x=   „    ,  x  distance  of  wing  a.c.      card  5. 

a.C. 

CTAIL  Ct,  Tail  chord,      card  6. 

2 

XKE  k   ,  radius  of  gyration  squared,      card  6. 

UE  LL^    o  relative  elevator  density,      card  6 

Bl  B,  ,  C,  .      card  6. 

1      bocr 

B2  B    ,  Ch         .      card  6. 

B2PRIME  B0,  1/2C,      .        card  6 . 

ZETA  X>>     ,  card  6. 

ZETADOT  }<   ,  card  6. 

AMETA  "*!*  ,  card  7 . 

AMETADT  >no  ,  card  7  . 

XME  m    (  mass  elevator,    card  7. 

e 

SE  Se,  area  elevator .    card  7. 
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Fortran      Program     LONGSJAV 

..JOE  BELL»LCNGST*E      1.5  MINUTES 

PROGRAM  LOKGSTAB 
C      PROGRAM  FCR  THE  CCMPUTATICN  OF  THE  RCC.TS  CF  THE  CHARACTER  IS<1  IC 
C      DYNAMIC  LONGITUDINAL  STABILITY  (STICK  FIXEC)  EQUATIONS 

CI  MEN  SIGN  KAPPA (10),CPR( 129),CPI( 129) .ROOTRt 1 28), ROOT  I ( 128) 

DIMENSION  CRR(129),CRI(  129  )  ,XR  {  M  )  ,  XI<  U  )  ,FX*  i*  )(,FXI  |U  ) ,  SR  (  3  )  JSf  (3  ) 

REAC  *6,  NOS,  MODE 
66  FORMAT  U3,I1    ) 

CO    777    J    *1,    NCS 

ICOUNT    «    C 
OCALL    INPUT1     <CLALFA,EPlA,C,TL,XCG,XAC,ST,SWiCLAtFATt.ETiATAIU$e»BLt 
1    RHO,G,fc,EPSALFA,THETA1,YH,CT,CMALFA.YI,V.€SLCBEfCSUfiL, KAPPA  ) 

C  COMPUTATION    OF    STABILITY    DERIVATIVES 

Q    *    ,5    «RHO    «    V«»2 

TCARAT    «    W    /<G«RHC«V*SW) 

THETAC  *  THETA1  /  57.2957795 

XU  ■  -CSLBD 

XW  «  .5«CSUBL«( l.C-t2.0«CLALFA/EPIA)) 

ZW  *  -.5»(CLALFA  ♦  CSUBD  ) 

APRIME  *  ,HU    ♦  2.CM  (XAC-XCO/C) 

BPRIME  *  -.58  -  CLALFA  »  APRIME  *  ( (XAC-XCG) /C) 

WINGMC  *  0.25  *((C/TL)»«2)«8PRIME 

TAILMC  *  -.50*  CLALFAT  «  (ST/SW)  •  ETATAIL 

BODYMC  "  -O.Cl«    <SB/SW)«(<BL/TL)««2       ) 

TOTALMQ  =  WINGMC  ♦  TAILMC  +  BOCYMQ 

YIB  =  YI/l(W/G)»TL»«2) 

WINGZC  *  -.25»CLALFA«(C/TL)«{.Ul4+(XAC-XCG)«2.0/C) 

TAILZC*  -.5*CLALFAT«(ST/SW)«ETATAIL 

Ul  *  */lG«RFO«SW»TL) 

BARMWCT  *  -.5»CLALFAT*ETATAIL«(ST/SW)«EPSALfA 

TOTALZQ  *  WINGZC  ♦  TAILZQ 

ZU  =  -C5UBL 

XTHETA  *  -.5»CSUBL 

ZTHETA  *  -.5*CSUBL«T/>NF(THETAC) 

AMU  *  -CT  »VH/TL 

AMW  *  .5  «  CMALFA  *  C/TL 
PRINT  5C5,ST,SW,  CLALFAT,ETATAIL,SB,BL,RHO,G,W«>EPSAtFA,THETAl,YH, 

1  CT.CMALFA.YI 

505OFORMAT  (1H1  /MCH         INPUT  PARAMETERS  // 

111H  TAIL  AREA=  F8.M.2X11H  WING  AREA=  F9.t4,2X,1UH  CL  ALFA  TML* 

2  F9.6,2X,1CH  ETA  TAIL*  F9.6,2X,1lH  BODY  AREA*  F9.2   / 

3  13H  BGCY  LENGTH*  FB.5t2X,  5H  RHO=  F1C.7.2X,  5H  GEE=  F8.5,2X, 
>4    8H  WEIGHT*  F12  .U,  1X,3HLBS  2X.18H  D  EPSILCN/D  ALFA=  F9.6   / 

5   21H  INITIAL  PITCH  ANGLE*  F9 .6 ,  1X.7HGEGREES  2X.22H  THRAJST  LINE  CI 

6STANCE*  F8.5,  1X,4FFEET  2X,16H  THRUST  CCEF  CT*  F6.6   / 

7    BH  CMALFA*  F9.6,2X,1CH  I  SUB  YY*  F  12. 3  ,  IX,  15HSLUG  FT  SCUAREC//) 

WRITE  OUTPUT  TAPE  5,  101 
1C1  FORMAT  (///  M7H  STABILITY   DERIVATIVES         // ) 

WRITE  OUTPUT  TAPE  5,  1C5  ,V, TCARAT, C 

105  FORMAT  (19H  FCRWARC  VELOCITY*  F12.3,1GH  T  CARAT«   F1C.6,  MH  C« 
1      F  1M.U      /         ) 

WRITE  OUTPUT  TAPE  5,  lC6t  CSUBL,  CSUBC 

106  FORMAT  (  9H  C  SUBL=    FK.6,  9H  C  SUBC=    F1C.6     /) 
kRITE  OUTPUT  TAPE  5  ,  1 C2,XU, XW.XQ, XTHETA , ZU.Zh ,TCTALZC 

1020FORMAT  (5H  XU*  F7.l4,6H   XW=  F7.5,6H   XC=  F7.U,1CH   XTHETA*  F7.5, 
1  6H   ZU*  F7.5,  6H   ZV«=  F7.M,12H   TOTAL  ZC=  F7.14    /) 
WRITE  OUTPUT  TAPE  5 , 1 C3, ZTHETA, AMU, AMW, BARMWCT,TCTALMC ,U  1  ,Y I E 
1030FORMAT  (  9H  ZTHETA*  F7.U,  6H   MU=  F8.6,  6H   Nfe*  F7.5,1UH   MEAR  W  C 
10T*  F7.M,  6H   MC»  F7.1*,  9H   MUONE*  F7.3.13H  I  SUB  BET**  F7.U   //) 
OCALL  PRERCOT  (XU, XK, >Q, XTHETA .ZU,ZW,TCTAL ZCl ZTHETA. AMU, AMW, BARMfcCT 
1, TOTALMCUl. YIB, TCARAT, Al,Bl, CI, Dl, El,  CRI,CRR     ) 
ROUTH  *  CRR(2)«  CRR(3)»  CRR(U)-  CRR(1)«  CR* ik )**2-CRR (2 )««2«CRR ( 5  ) 
WRITE  OUTPUT  TAPE  5,  108,  ROUTH 
108  FORMAT  {////35H  THE  VALUE  OF  ROUTHS  DESCRIMINANT  *   F10.2        ) 
IF  (RCUTH)  11,11, 12 

11  WRITE  OUTPUT  TAPE  5,  109 

1C9  FORMAT  <28H   THE  AIRCRAFT  IS  UNSTABLE      /?/   ) 
GO  TO  13 

12  WRITE  OUTPUT  TAPE  5,  110 

110  FORMAT  (26H   THE  AIRCRAFT  IS  STABLE        41/       ) 

13  CONTINUE 
"  N  ■  H 

CALL  FOLYMUL  (CRR ,CRI ,N, KAPPA, TCARAT) 

APPENDIX  II,  Fig.  2 
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DECISION   FCR  LONG  CR  LATERALt  STICK  FIXEC  {999  )f,  OR  STICK  FPEEC2C1) 
GO  TO  (777,201,  202  ,2CM ) ,MCDE 
2C1    CCNTINLE 

WRITE   CLTPUT   TAPE  5,301 
3010F0RMAT(////  80H     STICK  FREE  LONGITUDINAL  CYNAIUC  STABILITY  CCEFF 
1ICIEN1S  ////  ) 

WRITE  OLTPUT  TAPE  5f,  302,  KAPPA 
3C2  FORMAT  (  1CA8  ) 

CCALL  INPLT  2  ( CTA IL,XKE*UE,B 1, B2, B2PRIME ,ZETA ,2ETADCT, AMETM 
1      AMETADT    ) 
XKE  ■  SCRTF(XKE) 

IQ    =  TCTALZC 
AMQ  -  TCTALMQ 
OCALL  PRERCT2  (CTA  ILf. EPSALFA,  XKE,  Y  IB  ,U  1  ,U£,TL  ,  EHE2  ,82PR  IME,ZV, 
1  ZTHETA, ZQ,ZETA,ZETADCT,AMW,AMGtBARMWCT,AME<TA,ANETACT,  CRR.CRI    ) 
N  =  5 

CALL  FCLYMUL  (CRR ,CRI  ,N, KAPPA , TCARAT) 
777    CCNTINLE 
END  FILE  5 
ST0P2 
202    STOPU 
2C3    ST0P5 
20M   STOP* 
END 

CSUBROLTINE  INP0T1  (CLALFA , EPIA ,C ,TL, XCG , XAC ,ST ,SW .CLALFAT , ETATA  IL  . 
1  SB,.BL,RFO,G,W,EPSALFA,TFETAl,YH,CT,CMALFA,YI.V,CSUBC,CSUELtKAFf  A) 

DIMENSION  KAPPA(IC) .CPR(129),CPI( 129) ,ROOTR ( 1 28 ) ,RCCT I ( 128 ) 

CIMENSICN  CPR(  129  )  ,CR  I  (  1  29  )  ,XR  (  M  ) ,  XI  (  M  )  ,FXP  CM  )',FX  I  (M  ) ,  SR  (  3  )  ,  S  I'(3  ) 

READ  65*  KAPPA 
1CM  FORMAT  (  8F10.C  ) 

READ  10M,ST,SW.CLALFAT,ETATAIL.SB,BL,RH0,G 

READ  KM,  W,  EPSALFA, THETA 1 , YH ,CT,CMALFA  ,  Y  I  ,V 

READ  1CM,CSLBD,CLALFA,EPIA,CSUBL,C,TL,XCG,XAC 
65  FORMAT  (  1CA8) 

RETURN 

END 

SUBRCLTINE   PREROCT  (  XU,X  W,  XC,  XTHETA,  ZL  ,ZW,!ICTALZQ  ,ZTHETA,  AML't  AKW, 
ieARMWCT,TOTALMQ,Ul,YIB,TCARAT,A1.Bl,Cl,01,El,  CPI,  CRR) 

DIMENSION  KAPPA (10).CPR( 1  29  )  ,CP  I  (  129) , PCCTP ( 1 26 ) ,RCCTI( 128) 

01  MENS  I CN    CRR(129),CRI(129),XR{M),XI(M),FXR{M)tFXI(M),SR(3),SI(3) 

1  Al    *     l.C 

ZQMU    =    (1  .C    -MTCTALZC    /    V  1 )  ) 

2  El    =    -(XU    +    ZW)-t TCTALMC/YIB)-(ZQMU    «    EARNfcCT/Y18     ) 

3CC1    «    (XL*    ZW    -    XW«    ZU    )    ♦    (TCTALMC/YIE)«(XL*    ZW)    ♦  (    XL'«ZCMU   -    ZU    • 
KXQ/U1)    -    ZTHETA)*     (  BARMWDT/Y  IB  )    -(Ul «AMW/YIB  )» ( ZQMU    ♦    XC/L'l) 
.    MODI    =    -(TCTALMC/YIB)*    (XU»ZW    -    XW    •    ZU  )    *(-XTHETA«    ZU    ♦    ZTHETA    « 

1  XU    )«(BARMWDT/   Y  ie ) *(U1*AMW/ YIB )♦ (    XU«    ZCMi    -ZU«    XC/Ul    -ZrtHETA)* 

2  (U1»AML7YIB)»(-XTHETA   -XW»    ZCMU    ♦    ZW    «XQ/L1     ) 

50E1  =  (Ul«AMW/YIB)«(-XTHETA«ZU  ♦ZTHETA«XU)-(U 1*ANU/ Y IB ) ^-XTHETA* 
1ZW  +  ZTHETA*  XW  ) 
WRITE  OUTPUT  TAPE  5,1C1,  A1 ,B1 ,C1, D1.E1 
1C10FCRMAT  (/3H  A*  1X.Fl3.5tMH   B*  1X,F13.5.  MH   C*  1X,F13.5»  MH   C« 
1  1X,F13.5,MH   E«  1X,F13.5  //) 
CRR(1  )  -    Al 
CRR(2)  *  Bl 
CRR(3)  *  CI 
CRR(M)  *  D1 
CRR(5)  »  El 
,  •  CRKU  ■  0.0 
CRK2)  =  C.C 
CRK3)  '    C.C 
CRKM)  *  O.C 
CRK5)  =  O.C 
RETLRN 
END 

OSUBRCLTINE  INPUT2 (CTA IL *XKE ,UE,B1 ,B2, B2PRIM€, ZET A, ZETACCT , AM ETA, 
1      AMETADT    ) 

SUBROUTINE  BRINGS  IN  STICK  FREE  FACTORS 
CREAC  ltCTAIL,XKEtUEtBltB2tB2PRIME,ZETA,ZETACCT{tiAMETA,AMETACT,XME, 

1  FORMAT  (  8F1C.C  /MF1C.C) 
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WRITE  OUTPUT  TAPE  5,2 
2CF0RFAT(////1CX70H  INPUT    PARAFETEPS    STICK   FREE 

1  ) 

END 

SUBROITINE  PREROT2  (CTAI L  ,  EPS  ALFA, XKE , YIB .UI.Uf.TL.E 1 ,E2, E2PP  IFE, Z 
1W, ZTHETA,  ZQ,ZETA,  ZETACOT , AFW, AFC, BARMWCT ,  A>ET A, AFETACT , A, E         J 
;      SUBROITINE  CALCULATES  STABILITY  STICK  FREE  GUINT IC  COEFFICIENTS 
CIMENSICN         A(129),  E<129) 
UFACTCR  ~    2.0«  UE  / (U  1**2 ) 
UFACT  *  LFACTCR 
FACTOPK  =(XKE  /  TL)«*2 

XKFACT  -  FACTORK 
CTFACT   =  CTAIL  /  TL 
ZFACT    =  (l.C  -HZC/L1  )  ) 
EPSFACT   =  1,0-  EPSALFA 

WRITE  OUTPUT  TAPE  5.2C01,  UFACT,  XKFACT  .CTFACT  ,ZFACT,.EPSF  ACT 
2001CFORFAT  (M9H  FACTOR  CALCULATIONS  // 

1   9H  LFACTOR=   Fl  *4. 7 ,   9H  KFACTOR=   F1U.7,  1CH  CTFACTCR*   FU.7t 
210H  ZCFACTCR=   FT  i<  .7 ,   16H  EPSALFA  FACTCR=   F1U.7  //) 

;      A    CCEFFICIENT   CALCULATION 

A(1)=LFACT  «XKFACT  «  YIB 
:      B    CCEFFICIENT   CALCULATION 

0A(2)*LFACTOR»  XKFACT»(  BARFWCT  »(ZETACCT»CTFACT/U1  -ZFACT )-AMC-ZW* 
1  YIB  ♦  AFETAOT*  CTFACT  )  -  B2PRIFE  »  CTFACT  «  Yl E  /  Ul 
I      C    CCEFFICIENT   CALCULATION 

CA(3)=LFACT  "XKFACT  *(   ZW  #AFC  4  AFWM  ZETACOT  *CTFACT-  U1MZFACT) 

1  )  4  EARFWDT  »(ZETA  -  ZTHETA)+  Ul  *  AFETA  -  ANETADT  «  CTFACT  »ZW  ) 

2  4  E2FRIFE  «  CTFACT  /  Ul«(  AFC  *  ZW*YIE  4  ZFACT  *  BARFWCT  )  -  E2  • 

3  YIB  -Bl*(  AFETADT  «CTFACT/U1  4  ZETAOCT  *CTFACT/U1  •  (  YIE«EPSFACT 
U       4  BARFWCT  /Ul)) 

I      D    CCEFFICIENT   CALCULATION 

C  A  ( It )  -    UFACTOR  *  XKFACT   «(   Ul*AMWM  ZETA-ZTHETA  )-  t1«AFETA*ZW  ) 

1  -B2PPIFE  *  CTFACT  /Ul  »(   ZW  «AFC  -  Ll*AFW»   ZFACT  -  ZTFETA  * 

2  BARFWCT  )4B2*(AMC  4ZW  •  YIB  4  BARFWOT  *  ZFACT )~B1 *< AFETACT«CTFACT 

3  MEPSFACT  *  ZFACT  -  ZW/L1)4  ZETACOT  «  CTFACT/U1*( AFW  -  AFC  « 

4  EPSFACT  )4  ZETA  «(BARFWCT/Ul  4  EPSFACT  »  YIB  )  4  AFETA  ) 
:      E    CCEFFICIENT   CALCULATION 

GA<5)  =  B2PRIME  «  CTFACT  /Ul  *U1  *  AMW«  ZTHETA  4  B2«(  ZTFETA*EARFWD 
IT   -  ZW  »  AFG  4  Ul  *  AFW  *  ZFACT  )  4  E1*C  ZETA  *(EPSFACT  •  AFC  *■ 

2  AMW)  4  AFETA  «(  ZW  -Ul»EPSFACT  •  ZFACT  )-  AFETACT  «  CTFACT  • 

3  EPSFACT  *  ZTFETA  ) 

F    CCEFFICIENT   CALCULATION 
0A(6)  -    Ul«  ZTHETA*  (B2*  AFW  -  AFETA  •  B1«EPSFACT  ) 
CO   ICC   J=l,129 
100  B(J)  -    CO 

WRITE  CLTPUT  TAPE  5  ,  3,CT A IL, XKE, UE,B 1. E2 
3CF0RFAT  (///  1 3H  TAIL  CHORC  =  F7.3,1X,  8HK  SUB  E*  F9.5,1X,9FFL  SUE 
1E=  F9.5,1X,8HB  SUE  1=  F9.5,1X,  8HB  SUE  2=  F9.5      ) 
WRITE  OUTPUT  TAPE  5,4  ,B2PR  IFE , ZETA, ZETACOT, AFETA ,  AFETACT 
14CFCRFAK/1UH  B  SUB2  PRIFE=   F9.5,lX,  1  1H  Z  XUB  £TA=   F9.5.1X, 
115H  Z  SUE  ETA  COT=   F9.5,1X,11F  M  SYB  ETA=   F9.5,lX,UH  F  ISUE  ETA 
2  00T  =    F9.5  ///) 

CI  =  UFACT  «  XKFACT  »(  Ul  »  AFW  *  (ZETA  -ZTHETA)-  Ul«ANETA*  ZW) 
C2  *  -B2PRIFE  /Ul  *  CTFACT*(ZW  «AMQ-U1*AFW  *ZFACT  -ZTHETA,«EAPFWCT) 
C3  =  E2»(AFC  4ZW*YIB  4  BARFWCT*  ZFACT) 

004  =  -B1*(AFETA0T*CTFACT»<EPSFACT«ZFACT-ZW/UI )4ZETACCT/U1»CTFACT«( 
1  AMW  -  AFC*EPSFACT)  4ZETA*C  EARFWOT/Ul  4EPSFACT*YI B )  4AFETA) 

C  =  01402  4C3  4  DM 
A(»4)  *  C 
60FORFAT  (//    UH  Dl=  F15.8,  UH  C2=  F15.8,  kh    03*  F15.8.UH  D«4«  F15.8 
1  .  9H  D  TCT*l =  Fl 5.E,       //   ) 

El  =62PRI       ^TFACT  *<U1*AFW»ZTHETA) 

E2  =B2  •{  BARFWCT  -ZW*AMQ  4U1  »AFW«  ZFACT  ) 

0E3  =81  *CZEU*(EPSFACT«AMQ  -AFW )4AMETA*<ZW-U 1 *EPSFACT«ZFACT )- 
1   AMETACT»CTFACT»  EPSFACT  «  ZTHETA) 

E  =  El  4  E2  4E3 
7CF0RFAT  {//  UH    E 1=  F15.8,UH  E2*  F15.8.MH  E3=  F15.8,  9H  E  TCTAL* 
1  F15.E,  //) 

OCT  *UFACT  *XKFACT  «(ZW  "AFQ4AMW* ( ZETACCT*CTFACT-U1 «ZF ACT )4EARMWCT 
1    *(ZETA-ZTHETA) 4U1*AMETA-AFETADT  «CTFACT«  ZW  ) 

C2  »B2PRIFE/U1«CTFACT«<AFQ  4ZW*YIB  4ZFACT*  €ARFWCT) 
CC3  *  -B2  «YIB  -B1MAFETADT/U1  «CTFACT4  ZETACCT/U1*CTFACT«( YI E  * 
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1   EPSFACT  ♦  BARMWCT/U1)  ) 

C  =  CI  *  C2  ♦  C3 

WRITE  OUTPUT  TAPE  S,  e,  CI.  C2,  C3,  C 

WRITE  OUTPUT  TAPE  5,6 ,01 , C2  .D3  ,OU ,0 

WRITE  OUTPUT  TAPE  5,7.  E1,E2,E3,E 
8CF0RMAT  (///  UH  C1  =  F15.8,  UH  C2*  Fl5.e,  UH  C3  =  *15.e,9H  C  TOTAL* 
1   F15.8,        //  ) 

RETURN 

END 

SUBROUTINE    POLYMUL    <CRR,CRI ,N ,KAPPA,TCARAT ) 

CIMENSICN    KAPPA(1C).CPRM29),CPI(  129J,R0QTPJ12e],RCCTlM2e] 

CI  MENS  ION    CRRC129).CRni29),XR<uWxi{U),FXRN>.PXlU),SRT3)fcSIl3> 

IMAX    *    5C 

DEL    *    C.l 

NUM  *  5 
RATIO  =5.0 
ALTER  =1.CCCG01 
EP1  -.CCC000OO0OOCCCCG000O1 
EP2  =.CCC00C01 
EP3  =.0CCCGC1 
EPU  =.CCCCCC01 
SR(  1)*  -.5 
SI(1)=  CO 
SR(2)*  C.5 
SI(2)=  CO 
SR(3)-  CO 
SU3)*  CO 
MODE  =  1 
NNN  ■  N-fl 
7U9  IF    (CRP(D)   750,751,750 

751  CO  753  L=1,NNN 
M  =  L  +1 
CRR(L)  *  CRR(M) 

753    CONTINUE 

752  N  =  N-l 
GO  TO  71*9 

750    CONTINUE 

65  FORMAT  (1CA8) 

DO  5  1=1, 1C 

KON  =  EI- 

IF(K*PPMI)-KCN)2,5,2 

5  CONTINUE 
GO  TC  999 

2       PRINT    65, KAPPA 

66  FORM*T(2X,2HN=IU,12X,5HIMAX=I3,1CX,UHNUM=I3,llX,UI-DEL  =  E8.2* 
661    6X,6FRATIC  =  E8.2,  iiX,6HALTER=F1C8/ 

6  62    2X,UFEPl=E8.2,6X,UHEP2=E8.2,6XfUHEP3=E8.2,6X,UHEPU*E8.2/> 
MOCE=MCDE-H 

67  F0RM*T(6H   SRI*  E12.5.6H   SI1*  E12.5.6H   SR2*  E12.5.6H   512* 
671  E12.5,6H   SR3=  E12.5,6H   SI3*  E12.5) 

NP1=N*1 
71GF0RMAT  (U6HCTHE  COEFFICIENTS  OF  THE  GIVEN  POLYNOMIAL  ARE     /■// 
1   11H  REAL  FART     ) 

501  WRITE  OUTPUT  TAPE  5,71 

6810F0RMATC11CH       A  B  C  £  iE 

1  F  G  H  /) 

687   CONTINUE 

WRITE  OUTPUT  TAPE  5,  680  ,  (CRR ( J  ), J*1,NP1 ) 
680  FORMAT  { 1H   8F15.7   /  ) 

WRITE  OUTPUT  TAPE  5,  682 
682  FORMAT  (16H  IMAGINARY  PART   ) 

WRITE  OUTPUT  TAPE  5,  680,  (CRU  J  )  ,  J*l  ,  NP1 ) 

502  DO  6  1*1, NP1 
CPR(  I)=CRR(I) 

6  CPU  I)=CRI(I  ) 

CALL  C0MAG(CRR(1  )  ,CR  I  <1 ) ,C1 .KE ) 
DO  3C5  K*1,N 
GO  TC(9, 609), MODE 
53    FORMAT(ll-02XUHRCOT5X7HITERANT6X2(  10HREAL   PART  1CX1CHIMAG   PART1CX) 
1,1CHRC0T  IS  IN/2X6HNUMBERUX6HNUMBER7X2C9H  CF   *COT 1 IX) »2( 9H  CF   R( 
2Z)1lX).1lHA  RADIUS  CF) 
9     PRINT  53 
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6C9   1=1 

NPl=N+2-K 

NN=NP1-1 

IF (K+l-N) 13.12.1  1 

11  CALL  DIVD(-CRR<2),-CRI(2>fCRR<1>.CRH1>,XR41)VXI(l>,KE) 
K1  =  1 

K2=1 

GO  TC  16 

12  AR=CPRC1) 
AI=CRI  (1) 
BR*CPR(2) 

cr-cpr!') 

CI-CRlIll 

K  1  =  1 

K2=l 
MU  =  2 

GOTO   123 
1U   XR(1)=CEARR 
XI(1)=CEARI 
GO  TC  16 

13  DO  IS  J  =  1,3 
XR(J)=SR(J) 

15  X  I ( J  )  =  SI ( J ) 
Kl  =  1 

K2=3 

16  M1  =  l 
M2=l 
M3=1 
MU=1 

700  DO  717  L=K1,K2 

701  ZR=XP(L) 
7C2   ZI=XI(L) 

CALL  PCLYNCM(NP1  ,2R  ,  ZI.CRR  ,CRI  ,RR,R  I  ,KE  ) 

710  FXR(L)=RR 

711  FXI(L)=PI 

CALL  CCNAGCRR.RI .PKAG.KE) 
713   RAC=ALTER*<PMAG/C1)««(1.C/FLGATF(NN) ) 
71U   GO  TC  (715.19, 179), F3 

715  GO    TC(716,718),MCCE 

55      FORM/T(  2(16,    MX),      ME2C11,      5X,       E1C.il       ) 

716  PRINT    55,    K, I.ZR, Z I , RRf RI ,RAC 
718       POLYC=PCLYN 

POLYN=PPAG 
IF(P*AG)717,300,717 

717  1=1+1 
IF(K*1-N)17, 300,300 

17  GO  TC  (  18,20C),M1 

18  VAL=CEL«PCLYN 
DBARP=RAD 
DBARI  =  C0 

IS  K1=U 
K2=U 
Ml=2 
M3=l 

101  ABARR  =  XR(1)-XR(3  ) 

102  ABARI  =  XI  (1)-XI(3) 

103  BBARP=XR(2)-XR(3) 
1014   BBARI  =  XI(2)-XI(3) 

105  APIBR=XR( 1  )-XR(2  ) 

106  AMIBI=XI(1)-XI<2  ) 

107  CALL    MLLT(ABARR,AEARI,BBARR, BBARI ,TA . TE, K£  1 ) 
1C8      CALL    PLLT(TA,TB, AM  ER,  ANIBI  ,DENR,DENI  ,KE2) 

CALL    CCMAGtOENR.DENI.TA.KE) 
CALL    CCNAGIXR(3)  ,XI  (3).Tl*,KE) 
IF(TA-EP1»TU)  110,110,111 

110  CALL    CERIV(NPUXR(3),XH3),CRRtCRIfDR,0I,K50) 
CALL    CCNAG(RR,RI  ,TR,KE) 

CALL    CONAGCCR.CI  ,TC,KE) 
IF(TR-EP14»TC)192,171,  171 

111  DELAR=FXR(1)-FXR<3) 

112  0ELAI  =  FXK1)-FXI  C3) 

113  0ELBR*FXR(2)-FXR<3> 
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114  DELBI=FXI(2)-FXI  (3) 

115  CALL  MULT(e6ARR, 66 AR  I , DELAR , CELA  I  ,TA , TB, KE5 ) 

116  CALL  NULT(ABARR,  Ae  AR  I  ,CEL6R  ,CELB  I,  TC  .TO  ,K;E6) 

117  CALL  CIVD(TA-TC,Te-TC,0ENR,0ENI,AR,AI,K£7) 

118  CALL  NULT(ABAPR,  AE AR  I  ,TC  ,TO ,T  1  ,T2 , KEE  ) 

119  CALL  MlLT(eBARR, 6BARI,TA,TB,T3,TM,KES) 

12C  CALL  CIVD(T1-T3,T2-T4,0ENR,DENI,BR,BI,KE10) 
CR=FXR(3) 

ci'Fxim 

123  CALL  MULT(BR,BI,BR,ei,Tl,T2,KE11) 

12U  CALL  MULT(AR,AI,CR,CI,T3,T4,KE12) 

121  TA=T1-U.C«T3 

152  TB-T2-»l.C»T4 

CALL  CSCRT(TA,TB,TC,TD) 
147   T1=-6R+TC 
1 1*8   T2»-EI*TD 
149   T3=-ER-TC 
15C   T4=-EI-TD 

CALL    CCPAG(T1,T2,TA,KE14) 

CALL    C0KAG(T3.TUfTB.KEl5) 

153  IF(T/-TE)154,l68,16e 
1514       TA  =  TE 

155  T1*T2 

156  T2«T*4 

168  GO    TC    ( 157,159), M4 

157  IF(TA)161, 161,156 

158  CALL    COMAG(2.0*CR,2.0«Cl,T8,KE16) 
IF(TE-RA0«TA)159,l59,18C 

159  CALL    CIVD(2.0*CR,2.C«CI,T1fT2,DBARR, 0BARI,KE17) 
GO    TC    (  161,1*4)  ,M4 

161  XR(4)=XR(3)+CBARR 

162  XI(4)=XI(3)+D8ARI 
TR=AESF(XR{4)) 
TI»AESF(XI(4)  ) 
IF(TP)167,  167, 140 

140   IF(TI)167t167t169 

169  IF(TR-TI  )164,  167,163 

163  IF(TI-EP2«TR)165,167,167 

165  XI(4)»C.C 

GO  TC(5C3,504),MCDE 

56  F0RM*T(40HCITERANT  ALTERED  TO  BE  PURE  REAL  NLMCER.) 
5C3   PRINT  56 

5C4   GO  TC  167 

164  IF(TR-EP2*TI)166, 167,167 

166  XR(4)=0.0 

GO  TC(5C5,  167),MCCE 

57  FORMAT(45HCITERANT  ALTERED  TO  BE  PURE  IMAGINARY  NUMBER.) 
5C5   PRINT  57 

167  GO  TC  7CC 

180  CALL  DIVD(T1,T2, TA  , C.C,T 1 , T2, K20) 
CALL  CIVD(CR,CI, TB  , C.O,CR,C I, K2 1 ) 
CALL  MULT (CR, CI, RAC ,C. 0,CR,C I,K22) 
GO  TC(5C6,507)  ,MCCE 

58  FORMAT(87H  ITERANT  IS  OUTSIDE  CIRCLE  WHICH  BOUNDS  A  ROOT.   INTERP 
SeiOLATE  ITERANT  TO  ECGE  OF  CIRCLE.) 

5C6   PRINT  58 
507   GO  TC  159 

171  CALL  CCMAG(ABARR,A6ARI,T1,KR1) 

ErA(l)     STA(ITA)   ENA(3)     STA(ITB) 
IF  <T1-EP3«T4)  174,1  7*4,172 

172  CALL  CCNAG(BBARR,E6ARI,T2,KR2) 

ENA(2)  STAUTA)       ENAC3)  STAUTB) 

IF(T2-EP3»T4)175, 175,173 

173  CALL    CCMAG(AMIBR,AMI8I,T3,KR3) 

ENA(l)  STAUTA)       ENA(2)  STA(ITB) 

IF(T3-EP3«T4)174,174,111 

174  ENA(l)  STA(Kl)         STA(K2)         ENA(2)  STACM3)         SLJ(17E) 

175  EISAC2)  STA(KI)         STAIK2)         ENA(3)  STACM3) 
178       XRS(K1)=XR(K1  )«(1  .C+2.0«EP3) 

XIKK1)  =  XI(K1)«(1  „C+2.0«EP3) 
GO    TC(5C8,5C9)  ,MCDE 
6C         FORMAT! llHCITERANTS    XII, 6H   AND    XII, 41H    ARE   ffOC   CLOSC    TOGETHER      ALT 
1ER    ITERANT    X    IV, 1H.) 
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5C8   PRINT  6C,ITA,ITB,ITA 
5C9   GO  TC  7C0 
179   PCLYC=PNAG 

GO    TC    19 
2CO       IFCPCLYN-VAD201 ,2C1,210 
2C1       VAL=CEL«POLYN 
2C2      LIM=I+NUM 
203      M2=2 

CALL    CERIV(NP1,XR<M) ,XI < M  )  ,CRR,CRI ,DR, 0  I  ,K60  ) 

CALL    CCMAG(RR,RI  ,TR,KE) 

CALL    CONAGCDR.OI  ,TC,KE) 

IFiTP-EPU*TC)192.21C,21C 

210  OLT«RATI0«P0LYO/PCLYN 

211  IF(1.0-CLT)22C, 220,212 

212  CALL  ^L'LT(CEARR,CeARI,CLT,C.C,DBARR,DEARIfK2C) 
215   LIM=LIK+1 

GO  TC  (510,511),NCCE 
72  FORMAT ( 1 20H0PCLYNCM I AL  HAS  INCREASEO  IN  MAGNITUDE  TOO  MUCH  WITH  CU 
721RRENT  STEP.  THEREFORE  REDUCE  CURRENT  STEP,  ) 

510  PRINT  72 

511  GO  TC  161 

220  GO  TC  (221. 231), M2 

221  IF(I-INAX)250,250,222 

222  PRINT  62 

CO  68»*  J=1,NP1 

684  CRR(J)  =  10. C  «  CRR(J) 
ICOUNT  *  ICCUNT  ♦  1 

IF  {  ICCUNT  -  3  )  685.685,777 

685  fcRITE    OLTPUT    TAPE    5,686 

6860F0RMAT    (    55HCTHE    REVISED    COEFFICIENTS    CF    TEH  *HVEN   POLYNOMIAL    ARE 

1  ///) 
K    *     1 

GO    TO    687 

62  F0RMAT<69HCMAXIMUM  NUMBER  OF  ITERATICNS  POCFEB  WITHOUT  REOUCING 
621P(Z)  eY  CELTA.) 

513   GO  TC  777 

777     RETURN 

231   IF(I-LIM)250, 250,300 

250  DO  2fl  L  =  l,3 
XR(L  )=XR(L+1  ) 
XI(L)=XI(L*1) 
FXR(L)=FXR(L*1) 

251  FXI(L)«FXI(L+1) 
GO  TC  101 

192   GO  TC(51U, 300), MODE 
69   F0RMAT(19HCFIRST  DERIVATIVE  (  E17.9, E20.9, S5H)  INDICATES  THAT  ITE 
1RANT  IS  SUFFICIENTLY  CLOSE  TC  ROOT.) 
5H4   PRINT  69, DP, CI 
300   DO  3C2   J=2,NN 

CALL  MULT(ZR,ZI,CRRU-1),CRI<J-l)tTR»TI,K4C) 
CRR(  J)  =  TR  +  CRRU) 
302   CRK  J)  =  TI+CRI(J) 

63  F0RMATIS9H0THE  COEFFICIENTS  OF  THE  REDUCED  POLYNOMIAL  ARE  AS    FCLL 
631CWS  ) 

ROCTP(K)=ZR 
ROCTI(J<)=ZI 
GO  TC(3C3,3C5),MCCE 

303  PRINT  63 

68   FORMATdH  6E18.9) 

304  PRINT  68,(CRR(J), J=1,NN) 

PRINT  66,  (CRKJ  ),J=1,NN) 

305  CONTINUE 

GO  TC  (515,516). KCCE 

515  PRINT  75 

516  PRINT  EC 

8GCF0RI"AT(  1HC42X 1 4HT ABLE  OF  RCCTS//3X4HRCCT8X   1CHREAL   PART10X10HMA 
1G   PART  10X   6HPERI0C  10X  26HTIME  TO  CAMP  TO  HALF  AMPL 

2  /2X6HNUNEER7X2(9H  CF   R00T11X)  ) 
DO  3C6  1=1, N 

CALL  PCLYNCM(N+1  .RCOTRC I  )  . RCOT  I  (  I  ),CPR,CP'IVRR,R  I.  K»E  ) 
THALF  =  -.69314718C6/RC0TR(I)«  TCARAT 
PERIOC  »  6.2831853C72  /ABSFIROCTI ( I ) )  «  TCARAT 

306  WRITE  OUTPUT  TAPE  5,  84,  I,  RCCTRC  I  >, RCCTIU )  .PERICC , THALF 
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8U  FORMAT  (I6,UXtF16.6,3X1F16<6t7XtF10.llt8H  SECS     5X',F1C.U,5H  SECS) 
83   F0RMATU6,  MX,  ME2C.11) 

PRINT  75 
75   FORMAT(lHl) 
INAX  =  50 
NUM  =  5 
CEL  =  C.l 
RATIO  =  5.C 
ALTER  =1.C0C001 
EP1  =.CCCCCC0C00OC0CCCCCCC1 
EP2  =.CCC00001 
EP3  =.0000001 
£PU  -.0CCCCCC1 
SR(  1)*  -.5 
SI(  1)«  0.0 
SR(2)=  0.5 
SH2)=  C.C 
SR(3>=  CO 
SI  ( 2  )  *  0.0 
tfODE  =  1 
GO  TC  1 
1  RETURN 
999   CONTINUE 
END 

SUBRCUTINE    DER  IV( N,ZR, ZI , CR,CI ,DR,DI ,KER) 
OINENSICN  CR(129),CI(129) 

ENA(C)     STA(DR)    STA(CI)    STA(RR)    STA< RI ) EN A< 1 )STA (KER ) . 
DO  2  J=1,N 

CALL  FULT(ZR,ZI,RR,RI,TRR,TRI,K1) 
CALL  NULT(ZR,ZI, CR,CI  .TCR.TO  I, K2 > 

RSC(Kl)    AJPKL42)  RSC(K2)    AJP  (  1 )     ENA<2>  ST  Ai(KER  )SLJ  (2  )  . 

1  DR=TCR+RR 
DI-TCI*RI 
RR=TPR*CR( J) 

2  RI=TPI4CI(J) 

3  CONTINUE 
END 

SUBRCUTINE         PCLYNOMN.ZR.ZI.CR.CWRRfRI.KER) 
DIMEKSICN    CR(129),CI(129) 

ENA(O)  STA(RR)         STA(RI)         ENA  (1 )  STA(KER) 

DO    2    J=1,N 
CALL    NLLT(ZR,ZI,RR?RI,TR,TI.K1I 

RSC(Kl)         AJP(1)  ENA(2>  STA(KER)       SLJ(3) 

1  RR=TP+CR(J) 

2  RI*TI+CI(J) 

3  CONTINUE 
ENC 

SUBRCUTINE  CSCRT ( XR ,XI , YR, YI ) 
CCN(SQ2*1.i4lU2135624). 

LCA(XR)    LDQ(XI)    AJF21L+1)  LAC(XR)    CJP2(L+1)  LCC(XI) 
STA(A)     STQ(B)     CJPM1)    STC(C)     SLJM8)  +STA (P >SLJ ( S ). 

1  AJPK2)    LCA(B)     SLJU(8)   +FDVISC2)   STA(P) 
SSK(XI)    SLJU  +  2)   L*C(P)    ♦STA(C)     SLJ(5) 

2  TFS(B)     SLJ(3)     STC(S)     FCV{B)     STA(T)     ST A (R) SLJ (M) 

3  STA(S)     LDA(B)     FCV(A)     STA(T)     LCA(l.C)   STA <R ) SLJ ( U ) 

7  Y-SQPTF(X) 

8  SLJ(*)     STA(X)     SLJ(7) 

U      LCA(T)     FMU(T)     FAC(l.O)   SLJM8)   +F/IDCR)     FCV(2;C) 
SLJ*4(8)   4STA(T)     LDA(S)     SLJME)   «F*UT)     STA(P) 
LCA(XI)    FCV(2.C)   FCV(P)     STA(C) 

5  SSK(XR)    SLJ(6)     LDA(G)     LCC(P)     AJP2(L*1>  LAC(C) 
SSK(XI)    LDQ(P)     STA(YR)    STC(YI)    SliJ(L«-3) 

6  LCA(P)     LCQ(C)     STA(YR)    STQ(YI) 
ENO 

SUBRCUTINE   COHAG ( XR. XI »Z ,KER ) 

LCA(XR)         LDQ(XI)  AJP2(L*1)    LAC(XR)         CJP2(L*1)    LCC(XI) 

STQ(T)         *THS(T)  LLSCU8)         Q JP ( 2  )  STCCT)  FCV(T) 

STA(H)  FMU(H)  FAO(I.O)       STA(H) 

Y«SQPTF(H) 

-FfU(T)         ♦EXF7tlM1B)SLJ(L  +  2)       ENC  II )  SL\I43)  vERC(2) 

3  STQ(KER)       STA(Z) 

END 
SUBRCUTINE         OIVC (XR ,XI , YR, YI»ZR,ZI ,KER > 
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CAL 

2 

1   T*D 


ENC 

sue 

CAL 


ENC 
SUE 
CAL 
CAL 
PR* 
PI  = 
ENC 
SUB 


1A 
1 


ENC 
ENC 


L  PRC 
LCA(B 
ENA(2 
R«DR  + 
LCA(E 
LCA(P 
LCA(P 
STA(K 

RCUTI 
L  PRC 
LCAtB 
LCAiP 
LCA(P 
ENA(  1 
ENAC2 

RCLTI 
L  NCR 
L  NCR 
AP«CR 
AI«CR 

RCUTI 
SLJd 
LCA(  1 
QJP2( 
AJPK 
LCA(A 


C(XR,XI,YR\-YItBl,B2fPR,PIt0R1CI) 

2)         AJPI(l)         ENA(3)  SLJ(3) 

)  SLJ(2) 

CI«CI 

1)       -FDV(B2) 

R)         FOV(T) 

I)         FCV(T) 

ER) 


+  EXF7<ll4lB)SLJ(2)  STA(BI) 

-FNL(El)       +EXF7(1U1E)SLJ{2) 
-FNUCE1)       ♦EXF7(lUie)SLJI2) 


SIA(ZR) 
STA(2I)ENAC1) 


NE 


PULT(XRfXItYRtYItZRfZXffKER) 


0(XR,XI,  YP.YI, Bl,B2tPR,PI,OltC2) 
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I) 

) 

) 


-FMU(B]  ) 
-FNU(BI) 
-FMU(B1  J 
STA(KER) 
STA(KER) 


>EXF7(ll4lB)SLJ(l) 
SLJ(L*2> 


STMB1  ) 
STACZRJ 
STACZI) 


NE         PROC<XR,XIfYR,YI,BlfB2»PR,PI,CR,CI) 

K<XR,XI,81(,AR,AI) 

fMYR.YI, B2,CR,DI) 

-AI»0I 

♦AR«CI 


NE         NCRMA12A2,B1,S1,S2) 
)         ♦SEV7t70000E)  ZRC(C) 

A+1)    LCQ(Al)         CJP2(L+1)    LCC(Al) 
L*l)    LQC(A2)         LCL(1A*1)4THS(E) 
L*2)    STA(SI)         STA(El)         SLJU  +  5) 
1)         FDV(EI)         STACS1)         LDACA2) 


♦  ZPC<i»OOOE)ZRC<C) 
STL(E)  LCG<*2) 
SLJIL+2)       LCA(E) 

♦  ACCUA+2)     STA(El) 
FDV(BI)       *STA(S2) 


CC21 

RUN  1  CG  26    1 

93.14  375. C 

21CC0.0  .i<772 

.195  2.759 
F8U  LCNGITLCINAL 

93.4  375. C 

22C00.O  .14772 

.19  2.59 

RUN  3  C.G.  2k, 


93. U 

22COO.O 

k263 


375. C 
.11772 

2.92C 


'  22U  FPS   F8U 
2.393      l.C 
8.1        -.M37 
9.38U5     .902 

RUN2  STICK  FIXED 

2.392  .95 
5.6       -.»437 

9.38U5     .90 
V»139  KTS.23U  FPS 

2.393  1.0 
8.1                   -.M37 
9.38U5           .90 


1.29.63        BELL 
288.5 
C.195 
11.78 

M».209.V=23l4 
232. C 
.C511 
11.78 

3/3/63  LEVEU 
288.5 
0.195 
11.78 


RUBBERS 
52.8 
-.2937 
114.08 

CG2i»'*V« 
52.8 
-.380 
1*4. C8 

FLITE 
52.8 
-.380 

n.ce 


i. 


.002278 
;922CC.C 
29.625 
12VSL 
.CC2278 
96CCC.C 
29.375 

.CC227E 
96CCC.C 
29.375 


22. 

2214. 

29.! 

32." 
23*4 
29. ! 

22. 

23**. 

29.! 
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APPENDIX  III 
SYSTEM  MECHANIC  ATION  FOR  THE  ANALOG  COMPUTER 

Fig.  1  to  this  Appendix  displays  the  analog  schematic  for  the 
equations  of  motion  of  the  airframe.    Fig.  2  includes  the  schematic 
for  the  System  1  auto-throttle ,  aircraft  engine ,  and  stick  system 
dynamics  simulations.    Fig.  3  shows  the  System  2  auto -throttle 
schematic . 

Table  I  lists  the  potentiometer  settings  for  all  systems. 
Table  II  lists  the  scale  factors  used  in  computing  these  settings. 
Table  III  describes  the  calculations  used  in  arriving  at  these  values 
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Appen^iy  H£  TU'e  X    (s«n<w«/) 


TABLE  II 

SCALE  FACTORS  USED  IN  ANALOG  SIMULATION 

Scale  Factor  \  due 

b($.        =  .  005  radians/volt  =  .286  degree s/vol1 

^t>B-      =  „  005  radians/volt  -  .286  degrees/volt 

b/s*       =  .002  radians/volt  =  .115  degrees/volt 

^-Ook      =  .  005       tians/volt  =  .286  degrees/volt 

^Cu.      =  ,G02V    =.2%ofV    =  .278  Kts/volt 

o  o 

tfj^      =  .002  radians/volt  =  .115  degrees/volt 

°<r  =  200 1  :/  ,'oit 

C^TCi  =  40  lbs/radians  -  volt 

0(TC2  =  50  lbsAt  -  volt 

c*tc$  =  12    Lbs/radians -sec  -  volt 

^t>TCl  =  2:  ms  -  sec  -  volt 

n 

^jyrcz  -  3C  lbs/radians      sec    -volt 

C<0TC3  =  3C  ..    .  •■  ■>  -  sec^  -  volt 

&£&v  =  2,,:   k'         olt 

&tct  =  IOC       3/  '•olt  for  System  I ,         '     .    olt  for  System  2 

^repti6t  =  50  lbs/vo 

CK.Fi       =  2  lbs,  ■  c 

(XF3      =  =02  radians/volt 

o<DF3    =  .  1  lb,  sec  -  volt 

^fiiat  F     =  •'  '■"■  lbs/volt 


t*«     = 

.286  degrees/volt 

b<4r\y  = 

.02  g/volt 

Kw       = 

.01  ft/ft/sec  -  volt 
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/Appendry  HI       TaMc  IK 
Potent ow «^k    Self i ho,      Ca Icul&nons 


1. 

*«j_RCf  =  -.058 

/Al 

2. 

<Xtc»    K    =     .4 
o(tct 

Ifl 

3 

.           >«.  Ofe  R  =.  -.0(435 

/M 

1. 

Va.<v  |^  -.3433 

/Al 

5 

|M 

£ 

JM 

7. 

J  h 

8. 

**J*§-£-     -.,£41 

i  ai 

1. 

.       /i^t')o(i)^  r    ,<)9<9 

J  M 

10. 

^  *£.  S     -    .  ISO 

/M 

ii. 

^1-    l?a  =.  "2/39 

^6        ^t 

IM 

IE. 

»wA, ^  R  _  _/49 

,1/1 

13. 

1  M 

11 

1/1 

Rf.-.M 
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/^eJ/y    iH        I^Uf.    JZT       continue/ 

f5                          ^V^i^    =  -  3277  J  M 

1 7                    ^W  R  rf  -   .  loss'  I  n 

\8.                     ^m_  Ra    -   .7**  ll 

20.  '      2rsi  R     =.   *5oa  1/1 

(Xtct 

23.                      ML         =     .8^2  Rf*  '« 

21.  **  PC,  R    -     .6S45  J  A 

C0-  Oijci 

27                2lU-RG  =   .55'  sn     cf--.\/*f 

21              "          -L-Q<tc?  fi  =     .835  .5/1 

77    tf&TC? 


9  5 


30.  (^niR=.6oo  m 

31  ^HCt  =  .^2J  .5/1        C{--.l/*-f 

33,  >orci<!  -  .301  JM 

OC*v  fr 

35,  K2  *«  (I  =  ,388  I  Al 

2)6.  St/^fc     T^iMMfK  =,5"08 

0<F« 

38.  2L12L  K   =    ,/20  /0/c 

0(TCT 

4i.  A  K       -     .490         ,/ai 

4t.  ^*r»  R    =       £25" 


A  j? 

.490 

.270 

^Af 

«J 


e 
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Potent"!  ort\t^zy  C&lcuUf/OH 


is. 

48 


^i<^tR=  ,107 

tin  l 

fx*  J-<XihP=  -4+0 

fS        - X 


I  A) 
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Analog    ScAewat/c     /Airframe    E$uai'/ew    of    Moho^ 
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System  1     /Wfb-TVo/f/e     Ej^t/o*s  : 

AHZ  -  cha-hqe      of     norma,/     acce/erdT/0^    ^<. 


Fo 


W4S  4)4.  +  K 
Feel    S^ste^ 


Vdn^Ue  6a/.'n   =    3%£//W\     Stick: 
^  -  .0358   lb.  Set2  / if\ck 
b    -      /.?    /&.   /McA     st»&/<  /sec 
K  ~     20    lb.  /  l»c^     stick 

T  r    '/i?,r    sec 

-E   r     7.27    /isAdJ  /secz 
f 

tftdiy     JUL      F/q  ._2_      coh+mt/eg/ 

100 


To  fftjiM*. 


-loo 


AtJo3    MtcUtWU,    Sytcng     A^o-TUit/e 


f#U»tlmT.       /IT 


'  *4*-^[>**£] 


App€r\diy   1.1  E    .  Fi'q.-3 
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